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1.  Introduction 


This  is  the  final  project  report  on  the  research  program  entitled  "Elec¬ 
tromagnetic  Interactions  in  High-Speed  Integrated  Electronic  Circuits"  sup¬ 
ported  by  the  Office  of  Naval  Research  under  ONR  Contract  Number  N00014-86-K- 
0609.  It  reports  research  conducted  during  the  contract  period  of  June  1, 
1986  to  December  31,  1988. 


The  objective  of  this  research  was  to  investigate  electromagnetic  inter¬ 
actions  among  components  within  microstrip-based  circuits  fabricated  adjacent 
to  the  layered  surround  characteristic  of  a  high-speed  MMIC  environment.  EM 
fields  maintained  by  currents  immersed  in  layered  dielectric/ferrimagnetic 
media  are  exploited  to  study  the  complete  propagation-mode  spectrum  of  micro¬ 
strip  lines  and  their  interactions  with  the  coupled  microstrip  devices  which 
they  interconnect.  Seven  related  research  topics  were  investigated  in  connec¬ 
tion  with  this  program,  and  they  are  described  in  Section  4  of  this  report. 


The  seven  scientific  collaborators  who  contributed  to  the  research  pro¬ 
gram  are  identified  in  Section  2  of  this  report.  Section  3  details  biblio¬ 
graphical  information  on  the  eighteen  research  papers  which  were  precipitated 
by  the  investigations  supported  by  this  contract.  Finally,  Section  4  includes 
a  technical  description  of  the  project  and  results;  the  seven  sub-sections 
there  provide  details  of  research  on  the  major  topics  investigated  in  connec¬ 
tion  with  the  program.  Extensive  reference  is  made  to  sciectific  papers  from 
Section  3,  copies  of  which  are  included  in  Section  4. 
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2.  Scientific  Collaborators  on  Project 

The  following  individuals  have  participated  in,  and  contributed  to,  the 
research  program  during  the  contract  period. 


Faculty: 

1.  Dennis  P.  Nyquist,  Professor  of  Electrical  Engineering,  Principal 
Investigator. 

2.  Byron  C.  Drachman,  Professor  of  Mathematics,  Co-Principal 
Investigator. 

Doctoral-Level  Students: 

1.  Michael  J.  Cloud,  Graduate  Assistant,  Electrical  Engineering  Ph.D. 
Degree  awarded  in  December  1987. 

2.  Yi  Yuan,  Graduate  Assistant,  Electrical  Engineering  Ph.D.  Degree 
expected  in  June  1990. 

Master’s-Level  Students: 

1.  Eric  W.  Blumbergs,  Graduate  Assistant,  Electrical  Engineering  M.S. 
Degree  expected  in  June  1989. 

Undergraduate  Students: 

.1-.  Paul  F.  Havala,  Undergraduate  Assistant,  Electrical  Engineering  B.S. 
Degree  awarded  in  June  1987;  presently  completing  his  M.S.  Degree. 

2.  Darius  Adamozyk,  Undergraduate  Assistant,  Electrical  Engineering  B.S. 
Degree  awarded  in  June  1988. 
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3.  Contributed  Papers  Arising  from  Research  Program 

The  following  eighteen  scientific  papers,  reporting  results  of  this  re¬ 
search,  were  produced  during  the  Contract  period. 


[1]  B.  C.  Drachman,  D.  P.  Nyquist,  and  M.  J.  Cloud,  "Accurate  evaluation  of 
Sommerfeld  integrals  using  the  fast  Fourier  transform,"  Nat.  Radio 
Science  (USNCAJRSI)  Meeting,  University  of  Colorado,  Boulder,  CO,  digest 
p.  73.  Jan.  1987. 

[2]  J.  S.  Bagby  and  D.  P.  Nyquist,  "Radiative  and  surface-wave  losses  in 
microstrip  transmission  lines,"  Nat.  Radio  Science  (USNCAJRSI)  Meeting, 
University  of  Colorado,  Boulder,  CO,  digest  p.  226,  Jan.  1987. 

[3]  J.  S.  Bagby  and  D.  P.  Nyquist,  "Dyadic  Green’s  functions  for  integrated 
electronic  and  optical  circuits,"  IEEE  MTT-S  Trans.  35,  2,  206-210  (Feb. 
1987). 

[4]  M.  S.  Viola,  D.  P.  Nyquist,  and  B.  C.  Drachman,  "On  the  Sommerfeld  inte¬ 
gral  representation  of  the  electric  dyadic  Green’s  function  for  layered 
media, "  URSI  Radio  Science  Meeting,  Virginia  Tech,  Blacksburg,  VA,  di¬ 
gest  p.  278,  June  1987. 

[5]  M.  J.  Cloud  and  D.  P.  Nyquist,  ’’Complete  propagation-mode  spectrum  of 
microstrip  guiding  structures,"  URSI  Radio  Science  Meeting,  Virginia 
Tech,  Blacksburg,  VA,  digest  p.  120,  June  1987. 

[6]  E.  W.  Blumbergs,  D.  P.  Nyquist,  and  P.  F.  Havala,  "Integral  equation 
formulation  for  a  circular  patch  antenna  in  a  layered  environment,"  URSI 
Radio  Science  Meeting,  Virginia  Tech,  Blacksburg,  VA,  digest  p.  149, 

June  1987. 

[7]  M.  J.  Cloud,  "Electromagnetic  Interactions  in  Integrated  Electronic 
Circuits,"  Ph. D.  Dissertation,  Department  of  Electrical  Engineering, 
Michigan  State  University,  East  Lansing,  MI,  December  1987. 

[8]  D.  P.  Nyquist,  M.  S.  Viola,  M.  J.  Cloud,  and  M.  Havrilla,  "On 
Sommerfeld-integral  electric  field  kernels  for  microstrip-based 
circuits,"  Nat.  Radio  Science  (USNCAJRSI)  Meeting,  University  of 
Colorado,  Boulder,  CO,  digest  p.  112,  Jan.  1988. 

[91  M.  S.  Viola  and  D.  P.  Nyquist,  "An  observation  on  the  Sommerfeld  inte¬ 
gral  representation  of  the  electric  dyadic  Green’s  function  for  layered 
media."  IEEE  MTT-S  Trans.  36,  8,  1289-1292  (August  1988). 

[10]  C.  H.  Lee,  J.  S.  Bagby,  Y.  Yu^n  and  D.  P.  Nyquist,  "Entire  basis  MoM 
analysis  of  coupled  microstrip  transmission  lines,"  Nat.  Radio  Science 
(USNC/URSI)  Meeting,  University  of  Colorado,  Boulder,  CO,  digest  p.  179, 
Jan.  1989. 

[11]  Y.  Yuan  and  D.  P.  Nyquist,  "EFIE-based  perturbation  approximation  for 
coupled  microstrip  lines,"  Nat.  Radio  Science  (USNCAJRSI)  Meeting,  Uni¬ 
versity  of  Colorado,  Boulder,  CO,  digest  p.  180,  Jan.  1989. 
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[12]  B.  Drachman,  M.  Cloud,  and  D.  Nyquist,  "Accurate  evaluation  of 
Sommerfeld  integrals  using  the  fast  Fourier  transform, "  to  appear  in 
IEEE  AP-S  Trans. ,  April  1988. 

[13]  Y.  Yuan  and  D.  P.  Nyquist,  "Electric  dyadic  Green’s  function  for  layered 
media  with  dielectric/magnetic  contrast,"  to  be  presented  at  1989  URSI 
Nat.  Radio  Science  (USNC/URSI)  Meeting,  SanJose,  CA,  June  26-30,  1989. 

[14]  Y.  Yuan,  J.  Vezmar,  G.  King  and  D.  P.  Nyquist,  "Coupled  microstrip  tran¬ 
smission  lines:  full-wave  perturbation  theory  and  experimental  valida¬ 
tion,"  to  be  presented  at  1989  URSI  Radio  Science  (USNCAJRSI)  Meeting, 
SanJose,  CA,  June  26-30,  1989. 

[15]  D.  P.  Nyquist  "Dyadic  Green’s  functions  for  EM  fields  in  the  layered 
PC/IC  environment — representation,  properties,  and  applications,"  to  be 
presented  at  PIERS  (Progress  in  Electromagnetics  Research  Symposium), 
MIT,  Cambridge,  Massachusetts,  July  1989. 

[16]  D.  P.  Nyquist,  "Deduction  of  EM  phenomena  in  microstrip  circuits  from  an 
integral-operator  description  of  the  system, "  Proceedings  of  the  URSI 
Triennial  EM  Theory  Symposium,  Stockholm,  Sweden,  August  14-17,  1989. 

[17]  M.  J.  Cloud  and  D.  P.  Nyquist,  "A  note  on  the  mixed  potential  represen¬ 
tation  of  electric  fields  in  layererd  media,"  to  appear  in  IEEE  MTT-S 
Trans. ,  August  1989. 

[18]  C.  H.  Lee,  J.  S.  Bagby,  Y.  Yuan,  and  D.  P.  Nyquist,  "Entire-domain  basis 
MoM  analysis  of  coupled  microstrip  transmission  lines,"  submitted  to 
IEEE  MTT-S  Trans.  ,  February  1989. 
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4.  Technical  Description  of  Project  and  Results 

A  technical  description  of  the  seven  majoi  research  topics  investigated 
in  connection  with  the  project  is  included  in  the  corresponding  sub-sections 
to  follow.  All  bibliographic  references  are  to  scientific  papers  precipitated 
by  that  research,  as  listed  in  Section  3.  Reproductions  of  research  papers 
relevant  to  each  subsection  are  included  following  those  sections. 

4. 1  Electromagnetics  of  layered  dielectric/ferrimagnetic  media 

A  major  accomplishment  during  the  contract  period  was  the  development  of 
a  full-wave  theory  for  the  electromagnetics  of  layered  dielectric/ferrimag¬ 
netic  media.  A  dyadic  Green’s  function  for  the  electric  field  maintained  by 
electric  currents,  immersed  in  a  tri-layered  substrate/film/cover  surround 
having  general  linear,  isotropic  dielectric  and  magnetic  properties,  was  ob¬ 
tained.  In  the  special  case  where  the  substrate  layer  becomes  perfectly  con¬ 
ducting  (negative  infinite  imaginary  permittivity),  the  tri-layered  background 
becomes  that  appropriate  for  the  MMIC  environment.  This  formulation  and  its 
results  are  reported  in  papers  [3,4,9,13,15,17]. 

The  Hertz  potential  boundary  conditions  of  [3]  are  generalized  in  [13,17] 
to  allow  magnetic  permeability  contrast  between  the  various  isotropic,  homo¬ 
geneous  media  layers.  A  two-dimensional  spectral  analysis  subsequently  leads 
to  the  transform-domain  Hertz  potential  maintained  by  currents  immersed  in  the 
tri-layered  environment.  The  space-domain  electric  field  is  recovered  by  in¬ 
verse  Fourier  transformation,  which  leads  to  a  Sommerfeld-type  integral  repre¬ 
sentation  for  that  field.  Subsequent  complex  analysis  in  the  2-d  transform 
plane  provides  for  a  decomposition  of  the  field  into  discrete  surface-wave  and 
continuous  radiation-mode  spectral  components.  The  complete  propagation  mode 
spectrum  of  the  tri-layered  environment,  into  wnich  the  field  of  any  immersed 
currents  can  be  expanded,  is  consequently  identified.  The  surface-wave  con¬ 
tribution  arises  from  poles  of  implicated  spectral-domain  reflection  and  coup¬ 
ling  coefficients;  those  TE/TM  poles  of  the  tri-layered  media  are  located  by 
the  solutions  to  generalized  eigenvalue  equations  which  involve  the  dielectric 
and  magnetic  properties  of  each  layer  and  the  film  thickness.  A  specializa¬ 
tion  of  the  TM  eigenvalue  equation  leads  to  precisely  the  recently  reported  .M 
surface  waves  supported  by  a  grounded,  lossy  planar  ferrite  slab. 

The  source-point  singularity  in  the  Sommerfeld-integral  representation  of 
the  electric  Green’s  dyad  for  a  layered  IC  environment  was  studied  for  the 
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first  time  in  [4,9].  An  innate  principal  volume  and  associated  depolarizing 
dyad  to  properly  accomodate  that  singularity  were  identified.  A  new  mixed- 
potential  formulation  for  the  electric  field  maintained  by  currents  immersed 
in  the  cover  layer  of  the  tri-layered  IC  environment  was  exposed  in  [17]. 

That  formulation  avoids  the  source-point  singularity  which  arises  from  certain 
derivitives  of  the  corresponding  Hertz  potential,  resulting  in  Sommerfeld  in¬ 
tegral  representations  which  are  relatively  rapidly  convergent. 

The  electromagnetic  analysis  of  tri-layered  substrate/film/cover  media 
provides  the  basis  for  study  of  electromagnetic  interactions  in  MMIC’ s.  It 
provides,  through  Hertz  potential  and  electric  Green’s  dyads,  a  concise  and 
elegant  description  for  EM  excitation  of  such  an  integrated-circuit  environ¬ 
ment,  in  which  no  phenomena  have  been  neglected.  If  the  film  layer  is  ferri- 
magnetic  with  moderately  high  permittivity,  new  physical  phenomena  are  expect¬ 
ed  since  the  surface-wave  regime  may  couple  with  important  microstrip  propag¬ 
ation  modes  unless  the  film  is  very  thin.  Each  of  the  subsequent  research 
topics  exploits  this  EM  analysis  of  the  tri-layered  environment;  its  conse¬ 
quences  and  physical  implications  are  consequently  investigated  in  connection 
with  each  topic. 
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Dyadic  Green’s  Functions  for  Integrated  Electronic 
and  Optical  Circuits 

JONATHAN  S.  BAGBY.  number,  ieee,  and 
DENNIS  P.  NYQUIST,  member,  ieee 

Abstract  — Layered  stnictuies  play  an  important  role  in  both  integrated 
microwaTe  circuits  and  optical  integrated  circuits.  Accurate  prediction  of 
device  behavior  requires  evaluation  o(  fields  in  the  system.  An  increasingly 
used  mathematical  formulation  relies  on  integral  equations:  the  electric 
field  in  the  device  is  expressed  in  terms  of  the  device  current  integrated 
into  an  electric  Green’s  function.  Details  of  the  development  of  the 
specialized  Green’s  functions  used  by  various  researchers  have  not  ap¬ 
peared  in  the  literature.  We  present  the  development  of  general  dyadic 
electric  Green’s  functions  for  layered  structures;  this  dyadic  formulation 
allows  extension  of  previous  analyses  to  cases  where  currents  are  arbi¬ 
trarily  directed.  The  electric- field  Green’s  dyads  are  found  in  terms  of 
associated  Hertzian  potential  Green’s  dyads,  developed  via  SommerfeM’s 
classic  method.  Incidently,  boundary  conditions  lor  electric  Hertzian 
potential  are  utilized;  these  boundary  conditions,  which  have  been  a  source 
of  confusion  in  the  research  community,  are  developed  in  full  generality. 
The  dyadic  forms  derived  herein  are  reducible  in  special  cases  to  the 
Green’s  functions  used  by  other  workers. 

I.  Introduction 

Layered  dielectric  structures,  such  as  those  depicted  in  Fig.  1, 
play  an  important  role  in  both  integrated  electronic  circuits  and 
integrated  optical  circuits.  In  integrated  electronics,  conducting 
“devic  .s”  are  affixed  to  a  dielectric  film  layer  which  is  deposited 
over  a  conducting  ground  plane.  For  integrated  optical  circuits,  a 
dielectric  waveguiding  region  is  typically  placed  on  top  of  a 
dielectric  film  layer;  the  film  layer  is.  in  turn,  deposited  on  a 

Manuscripl  received  July  14,  1986;  revised  September  6.  1986. 

J.  S  Bagby  is  with  the  Department  o(  Eiectrical  Engineering.  Uruversity  of 
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D  P.  Nyquist  is  with  the  Department  of  Electrical  Engineering  and  Systems 
Science,  Michigan  Sute  University,  East  Lansing,  MI  48824. 
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(a)  (b) 

Fig.  I.  Typical  structures  for  (a)  integrated  electronic  circuits  and  (b)  in¬ 
tegrated  optical  circuits. 
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Fig.  2.  Multilayered  dielectric  background  sinicture  with  source  in  <  th  iaver 


n';  the  primary  potenlial  propagates  directly  from  the  source  to 
a  field  point  in  the  ith  layer,  whereas  the  scattered  potenual 
arrives  at  a  field  point  after  being  scattered  (reflected  or  trans¬ 
mitted)  from  interfaces  between  adjacent  layers.  These  potentials 
satisfy  the  following  Helmholtz  equations: 


substrate.  In  both  cases,  accurate  prediction  of  device  behavior 
requires  evaluation  of  the  fields  in  the  system. 

Exact  formulations  for  the  electric  field  in  terms  of  differential 
equations  [1],  [2]  are  rendered  ineffective  for  these  more  com¬ 
plicated  geometries  by  the  inseparability  of  boundary  conditions 
at  the  device  boundaries,  and  approximations  based  on  such 
formulations  have  been  found  to  neglect  important  radiative 
phenomena  [3],  (4|.  Alternative  formulations  have  been  gaining 
favor;  these  formulations  rely  on  integral  equations  to  describe 
the  behavior  of  the  system  I4]-(8].  The  electric  field  in  the  system 
is  expressed  in  terms  of  a  current  integrated  into  an  electric 
Green's  function;  boundary  conditions  are  incorporated  in  their 
full  generality  in  the  Green's  function.  However,  details  of  the 
development  of  the  Green's  functions  used  by  different 
researchers  have  not  appeared  in  the  literature. 

We  present  a  general  dyadic  formulation  for  the  electric  Green's 
functions  for  such  layered  structures.  The  dyadic  forms  derived 
here  allow  extension  of  previous  analyses  to  cases  where  normally 
directed  currents  exist,  such  as  in  thick  devices.  The  electric-field 
Green's  dyads  are  found  in  terms  of  associated  Hertzian  poten¬ 
tial  Green's  dyads;  these,  in  turn,  are  developed  by  an  extension 
of  Sommerfeld's  classic  method  [9],  with  components  expressed 
as  2-D  spatial  frequency  integrals  of  the  Sommerfeld  type.  Inci- 
dently,  boundary  conditions  for  electric  Hertzian  potential  are 
utilized;  an  adaptation  of  Sommerfeld's  [9]  development  of  these 
boundary  conditions,  which  have  been  a  source  of  confusion 
(compare,  e.g.,  [10]  and  [11]),  is  given  in  the  Appendix  for  the 
reader's  convenience.  Application  of  the  Green's  dyads  to  the 
analysis  of  integrated  microstrip  circuits  and  integrated  optical 
waveguides  is  discussed. 

II.  Prototype  Potential  Solutions 

In  this  section,  we  develop  prot.  pc  solutions  for  the  Hertzian 
potential  excited  in  layered  background  structures  by  impressed 
currents.  These  solutions  will  later  be  specialized  to  the  cases  in 
Fig.  I  by  application  of  boundary  conditions.  The  resultant 
expressions  for  the  potential  in  terms  of  the  impressed  current 
allow  identification  of  Hertzian  potential  dyadic  Green's  func¬ 
tions,  from  which  associated  electric- field  Green's  dyads  can  be 
found. 

Consider  the  situation  depicted  in  Fig.  2.  An  impressed  current 
J  (or  an  impressed  polarization  P-J/ju)  radiates  in  the  ith 
layer  of  an  infinite,  multilayered  dielectric  structure,  generating 
electric  Hertzian  potential  in  each  layer.  The  total  potential  in 
each  layer  is  the  sum  of  a  primary  part  fl and  a  scattered  part. 


(V^  +  A:/)| 


/  -  ' 
\  0.  all  /. 


(1) 


The  solution  for  the  primary  potential  in  terms  of  the  impressed 
current  is 


nf(7)~jG^(T\7') 


At-) 

ju(, 


dV 


(2) 


where  C'(r|r’) -exp(-  Jk,R)/dirR  is  the  free-space  Green  s 
function. 

We  solve  (1)  for  the  scattered  potential  by  Fourier  transforma¬ 
tion  on  spatial  variables  tangratial  to  the  layer  interfaces  Define 
the  transform  pair  F(r)  ** /(X,  y),  where  X  is  the  2-D  transform 
variable  X  —  (Jc  +  fi.  Then  we  have 


-0 


(3) 


where  pj  -  X^  -  k}  kj.  Equation  (3)  has  solutions 

ii7(X,>-) -W'(X)exp(  ?/»,>-)  (4) 

where  the  coefficient  1F(X)  is  determined  by  application  of 
boundary  conditions;  (2)  can  also  be  transformed  to  facilitate 
application  of  boundary  conditions: 


•'  JU(, 


(5) 


Here,  J{r)  «-►  j(\,  y)  and  g'(X;  y,  y')  -  exp(  -  p,\y  -  .v'D/2;), 
We  will  use  (4)  and  (5)  as  prototype  solutions  for  the  trans¬ 
formed  Hertzian  piotential  in  the  structures  of  Fig.  1.  Since  the 
coefficients  W(X)  will  depend  on  the  transformed  impressed 
current  through  (S),  inverse  transformation  of  these  expressions 
will  allow  identification  of  the  Hertzian  potential  Green's  dyads 
for  these  structu.es,  with  the  total  potential  given  by 

n,(r)-/C(r|r)-- - dV’.  (6) 

■'  /ur. 


Once  the  Hertzian  potential  has  been  determined,  the  corre¬ 
sponding  electric  field  can  be  found  as  £;“(VV  +kf)n|,or 

-  (VV- +  X/)jC(r|r')  - - dV 

r-  A^) 

- dV'  (7) 

J  jut, 
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Fig.  3.  Detail  of  poieniiai  componenis  in  (a)  integrated  electronic  circuit  and 
(bl  integrated  optical  circuit. 


where  the  electric-Held  Green’s  dyad  C(r|r')  is  defined  by 

G'(r\r')  ~Py. {vV  +kf)G(r\r')+LS(7-P).  (8) 

P.V.  indicates  that  spatial  integr^ion  must  be  performed  in  the 
principal  value  (P.V.)  sense,  and  L  is  an  appropriate  depolarizing 
dyad  [12]. 

III.  Microstrip  and  Optical  Circuit  Green’s  Dyads 

We  now  apply  boundary  conditions  to  the  prototype  Hertzian 
potential  solutions  above  to  find  spedHc  Hertzian  potential 
Green’s  dyads  for  the  cases  in  Fig.  1.  For  the  sake  of  brevity, 
some  purely  algebraic  steps  are  omitted. 

A.  Mtcrostrip  Green's  Dyad 

First,  consider  the  microstrip  integrated  circuit  background 
structure  of  Fig.  3(a).  The  transformed  potentials  in  the  cover 
medium  and  film  layer  decompose  into  principal  and  reflected 
parts  as  follows: 

+  +  (9) 

where  c  and  /  refer  to  the  cover  medium  and  film  layer, 
respectively.  The  principal  potential  is  given  by  (5)  with  / 
replaced  by  c;  the  scattered  potentials  arc  written  in  terms  of  the 
prototype  solution  (4)  as 

w; -^exp(-;>.y)  e/‘ -  Bcxp(  -  Ccxp(  - 

(10) 

We  apply  boundary  conditions  to  evaluate  the  coefficient  /4(\); 
then  (9)  can  be  inverse-transformed  to  yield  a  general  expression 
for  the  potential  in  the  cover  medium  in  terms  of  the  impressed 
current. 

Application  of  tangential  Hertzian  potential  boundary  condi¬ 
tions  (see  the  Appendix)  at  >>  -  0  and  y—-t  results  in  an 
expression  for  the  tangential  components  of  A  in  terms  of  a 
tangential  reflection  coefficient,  R,(X) 

p,  -  p,  coth  p,t 

^  (11) 

p^  +  Pf  coth  pft 


where  the  source  term  y  is  given  by 

K‘f— - --zr-dy',  a~x,y,z.  (12) 

J  2p, 

Application  of  normal  boundary  conditions  at  y  -  0  and  ^  = 
-  t  provides  an  expression  for  the  normal  component  of  A  in 
terms  of  a  normal  reflection  coefficient  R„(\)  and  a  coupling 
coefficient  C(X) 


+C[7fl'.  +  yfK]. 

K/. 

_ 2(K-\)p, _ 

(  A  P/  coth p/i)(  Kp^  +  Pf  tanh p^i) 


Kp^  -  Pf  tanh  p/t 
Kp^  +  Pf  Pfi 


Having  completely  specified  A,  we  are  now  in  a  position  to 
solve  for  the  total  potential  in  the  cover  medium.  Inverse-trans¬ 
forming  (9)  after  substitution  of  (S)  and  (10)-(13)  gives  the  result 
in  dyadic  notation 

nAr)-fG(r\F')^-^dy  (i4) 

The  Heruian  potential  Green’s  dyad  in  (14)  is  found  to  decom¬ 
pose  into  a  principal  and  a  reflected  part:  C(r|r')  -  IG'’{r\r') 
-t-  G''(r\7').  with 

G'(  717' )  ~  xG.'x  +  y|  ^  C.'x  +  G^y  +  ^  C.'f  j  +  zG/z  (15) 

where  the  scalar  components  of  the  reflected  Green's  dyad  are 
given  in  terms  of  2-D  inverse  transform  integrals 

C;(r|r')\  ^  . . . 

G:(7i7')\ -  ff  (RJ\)\ - ^ - d-x  (16) 


y:(r]r')j-fj  (/?„(X) 
:i:(7\7')l  '’'lc(X) 


2(2ir)>. 


B.  Integrated  Circuit  Green's  Dyad 

Now  consider  the  optical  integrated  circuit  background  struc¬ 
ture  of  Fig.  3(b).  The  potentials  in  the  cover  medium,  film  layer, 
and  substrate  are  in  this  case 

€  =  +  =  (17) 

where  c.  /.  and  s  refer  to  the  cover,  film,  and  substrate,  respec¬ 
tively.  The  principle  potential  iff  is  as  in  (5)  with  i  replaced  by  c. 
The  scattered  potentials  are  written  in  terms  of  the  prototype 
solutions  (4) 

i^'  =  ;4exp(-p,>  )  if/'  •  BeKp(pfy) 

if/'  ~Cexpi-  pfv)  if;  •  Dexp(p,y).  (18) 

Once  A  is  found  via  application  of  boundary  conditions,  (17)  can 
be  inverse- transformed  to  give  a  general  expression  for  the  poten¬ 
tial  in  the  cover  medium. 

Application  of  tangential  boundarv  conditions  at  y  0  and 
y  “  -  I  once  again  gives  the  tangential  components  of  A  in  terms 
of  a  tangential  reflection  coefficient  R,{X) 


A  - 

A  +  A^: 


p,  +  pfUnhpft 


P,+P,l^i  '  Pf  +  P>^^P/i 

(19) 

where  the  source  term  A  is  given  by  (12). 

Application  of  normal  boundary  conditions  at  v  =  0  and  -  r 
yields  an  expression  for  the  normal  component  of  A  in  terms  of  a 
normal  reflection  coefficient  R„(X)  and  a  coupling  coefficient 
C(X) 


+C(y«r  +yfK]. 


7^1  A  -  A^< 
^1 A  +  A  ^4 


Pt  ^2  Pt  A  ‘ 
K-.  Pf  -t-  p,  tanh  p,  t 

(1^  A,)(l^  A,,) 

2f(,P, 


^i=VA, 


.  K~Cf/c,.  (13) 


I  l'^:P/  +  P.  j 

We  can  now  solve  for  the  Hertzian  potential  in  the  cover  by 
inverse-transforming  (17)  The  result  is  identical  in  form  to  (14). 
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The  Hertzian  potential  dyad  again  decomposes  into  a  principle 
and  a  reflected  part;  G(r\r')  -  /C'(Flr')+  C'(rlr'),  with  G'  as 
in  (IS)  and  (16)  and  R„  R,,  and  C  as  in  (19)  and  (20). 

IV.  Appucations 

In  this  section,  we  review  some  applications  of  the  electric 
Green’s  dyads  found  above  in  the  analysis  of  waveguiding  sys¬ 
tems  and  compare  them  with  specialized  versions  used  by  other 
workers. 


Pig.  4.  Hertziao  potential  at  a  dielectric  interface. 


A.  Microstrip  Integral  Equation 

Consider  the  microstrip  device  illustrated  iii^Fig.  1(a).  If  excita¬ 
tion  is  provided  by  an  incident  electric  field  a  surface  current 
K  is  induced  on  the  surface  of  the  microstrip  device.  This 
induced  current  in  turn  supports  a  scattered  field  E‘  in  the 
system.  The  tota^  tangential  electric  field  at  the  strip  surface  must 
vanish,  or  —t  E‘~  ?•£'.  The  scattered  field  may  be  written  in 
terms  of  the  electric  dyadic  Green's  function  derived  in  Section 
Ill-A  as 

K{  y* ) 

E'{r)~(G*{r\r')-^^dS'.  (21) 

J 

Thus,  we  obtain  an  integral  equation  for  the  induced  surface 
current  on  the  device 

r/G'(r|r')-^dS'--r£'(r).  (22) 

By  assuming  the  microstrip  device  to  be  infinitely  thin,  we 
recover  a  specialized  version  of  this  equation  that  has  been 
successfully  utilized  by  several  workers  in  the  analysis  of  thin 
microstrip  structures  [5]-(8].  This  more  general  version  provides 
the  basis  for  continuing  research  on  surface  wave  and  radiative 
loss  and  coupling  in  thick  microstrip  devices. 

B.  Dielectric  Waveguide  Integral  Equation 

Consider  the  integrated  dielectric  waveguide  of  Fig.  1(b).  If  the 
system  is  excited  by^  an  impressed  field  £',  then  the  total  field  in 
the  system  is  £—£'-(-£’,  where  £'  is  the  field  scattered  by  the 
waveguide  core.  The  waveguide  core  can  be  replaced  by  an 
induced  distribution  of  equivalent  polarization  ^~(c(r)- 
<,.)£( r),  and  the  scattered  field  can  be  written  in  terms  of  this 
equivalent  polarization  as 

£’(r) (23) 

Then,  setting  £  -  £'  •  £',  wc  have 

j  "7  —  r\r' )  ■  E(r')  dW  ~  E-(r)  (24) 

which  is  an  integral  equation  for  the  unlinown  total  field  in  the 
waveguiding  system.  This  equation  has  been  utilized  successfully 
in  the  analysis  of  propagation  and  field  distribution  in  integrated 
optical  waveguides  [4],  and  is  the  basis  for  continuing  research  on 
quantification  of  surface  wave  and  radiative  loss  in  integrated 
optics 


Appendix 

Hertzian  Potential  Boundary  Conditions 

Here,  we  present  an  adaptation  of  Sommerfeld’s  (9)  develop¬ 
ment  of  Hertzian  potential  boundary  conditions  to  assist  the 
reader  Consider  the  situation  in  Fig.  4.  Electric-type  Hertzian 


potential  exists  in  regions  1  and  2,  separated  by  a  locally  planar 
boundary.  The  fields  in  both  regions  are  found  in  terms  of  the 
potential^in  the  region,  it,  or  as  £- (A:^ -t- V  V  )n,  H- 
Ju€V  X  n.  Let  the  potential  in  each  region  at  a  point  at  the 
interface  have  components  11  -  rll,  -1-  all,,  with  r  tangent  and 
n  normal  to  the  interface,  and  b  —  txh.  The  fields  in  each 
region  are  given  by 

£,-/c^n,  +  -^(vn) 

fs--^(^  n)  £,-A:^n.  +  -^(v-n) 

d  id  3  \ 

w, -><->€— n,.  (25) 

Enfordng  continuity  of  tangential  £  at  the  interface  gives  V  -  111 
-V-flj  and 

n„-Arn,j  (26) 


with  K  “  (2/<i-  In  a  similar  fashion,  continuity  of  tangential  H 
yields 


5n„ 

Finally,  V  •  H,  «  v  •  flj  and  (26)  give 


dn 


d  ,  an„ 

^(n,.-n„2)-(A:-‘-i)^ 


(27) 


(28) 


Equations  (26)-(28)  are  thus  the  appropriate  boundary  condi¬ 
tions  for  electric  Hertzian  potential  at  a  material  interface. 

If  medium  2  is  a  perfect  conductor,  the  vanishing  of  tangential 
£  at  the  interface  yields  the  following  special  case: 

3n„, 

n,.-0  -^-0.  (29) 
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The  electric  dyadic  Green's  function  for  layered  media  is 
discussed.  It  is  well  known  that  for  the  free-space  electric  dyadic 
Green's  function  Gq,  evaluation  of  the  electric  field  at  observation 
points  withirj  the  source  region  requires  specification  of  a  depola¬ 
rizing  dyadX.  The  dyad  T  is  dependent  on  the  geometo  of  the 
"principal  volume"  which  excludes  the  singularity  of  Special 
considerations  are  invoked  for  the  layered  background  media  which 
are  appropriate  for  the  electromagnetics  of  integrated  electronics. 

The  relevant  equation  which  uniquely  defines  the  electric 
field  maintained  by  arbitrarily  located  electric  currents  in  a 
layered  environment  is  established.  Beginning  with 


E(r)»(k2  +  77- ) 


+  6'^)-J(?')dV' 


(1) 


and  using  the  ^tmierfe Id-integral  representation  for  the  dyads 
("principal"^  has  the  pertinent  singularity,  "reflected"^  is 
well  behaved)  appearing  in  the  integrand  of  (1),  it  is  shown  that 
the  electric  field  may  be  expressed  as 

E(r)=-ja)PO  Mm  f  G®(r;r '  )*J(r' )dV'  -  L*J(r)/jwe  (2) 
where  and  T  are  quantified. 

In  conclusion,  it  is  shown  that  for  the  layered-media  electric 
dyadic  Green's  function,  the  Sommerfe Id- integral  representation  is 
appropriate  and  leads  to  an  innate  option  for  the  principal  volume. 
The  correction  term  in  (2)  is  precisely  the  depolarizing  dyad 
corresponding  to  this  preferred  choice  of  the  excluding  region. 
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An  Obsenratioa  on  the  Sommerfeld-Integral 
Representation  of  the  Dectric  Dyadic 
Green’s  Function  for  Layered  Media 

MARK  S.  VIOLA  and  DENNIS  P.  NYQUIST,  member,  ieee 

Abstrmet  —The  electric  dyadic  Green’s  function  for  layered  diekctrics  is 
dbcttssedMt  is  i*eO  kaonrn  that  for  the  free  space  elec^  dyadie  Green’s 
function  eraluation  of  the  electric  field  at  observation  points  within  the 
source  region  requires  spedficaHon  of  a  “principal  vobune’’  along  with  the 
corresponding  depolarizing  dyad  L.  Special  considcralions  are  inroked  for 
layered  background  nMdia  which  ere  appropriate  for  the  dectronugnetics 
of  integrated  electronics.  It  is  shown  that  use  of  the  SounaerfeU-integral 
representation  of  the  electiic  dyadic  Green’s  function  leads  to  an  innate 
choice  for  the  depolarizing  dyad.  A  corresponding  prindpei  vohune  b 
subsequently  identified;  it  b  demonstrated  that  there  exbb  an  eltemative 
choice  for  thb  excluding  region  which  leads  to  the  same  depolarizing  dyad. 

I.  iNTRODUenON 

There  is  an  increasing  interest  in  the  study  of  optical  and 
electronic  circuits  inunersed  in  a  layered  dielectric  surround. 
Conventional  differential-operator  formulations  for  the  fields 
within  these  circuit  devices  are  rendered  ineffective  due  to  the 
inseparability  of  the  applicable  boundary  conditions  for  struc¬ 
tures  having  practical  shapes.  An  integral-operator  formulation, 
based  on  the  identification  of  equivalent  volume  polarization 
currents,  circumvents  this  difficulty.  Construction  of  the  integral 
operator  requires  knowledge  of  the  Green’s  function  for  the 
layered  surround. 

^A  general  development  of  the  Hertzian  potential  Green’s  dyad 
G  for  layered  dielectrics  has  been  discussed  by  Bagby  and 
Nyquist  (!].  Based  on  the  classical  method  of  Sommerfeld  [2],  the 
Hertzian  potential  dyadic  Green’s  function  was  shown  to  have 
scalar  components  represented  by  2-D  spectral  integrals.  As 
asserted  by  Yaghjian  [3],  the  singularity  of  G  is  seen  tojuise  from 
that  part  of  the  dyad  which  is  the  Green’s  function  G'  (“prin¬ 
cipal”  Green’s  dyad)  for  the  unbounded-space  problem.  In  Sec¬ 
tion  II,  the  development  in  [1]  is  alter^  slightly  so  that  identifi¬ 
cation  of  a  natural  depolarizing  dyad  L,  relevant  to  the  Green’s 
dyad  G'  for  the  electric  Geld,  may  subsequently  be  made  in 
Section  III. 

Finally,  the  electric  field  is  expressed  in  the^tandard  form  as  a 
volume  integration  of  the  scalar  product  of  G'  with  the  electric 
current  source  J.  The  volume  of  integration  extends  over  the 
support  of  the  current  density  but  excludes  the  singularity  point 
of  G'.  The  excluding  region  is  identified  as  the  “principal 
volume”  which  corresponds  to  the  preferred  choice  of  the  de¬ 
polarizing  dyad  as  tabulated  in  (3].  An  alternative  excluding 
region,  which  is  shown  to  be  equivalent,  is  suggested  to  be  useful 
in  practice  due  to  its  simple  form. 

II.  Hertzian  Potential  Green’s  Dyad 

In  this  section,  the  Hertzian  potential  dyadic  Green’s  function 
is  developed  for  the  trilayered  structure  depicted  in  Fig.  1.  A  film 
layer  of  thickness  /  and  refractive  index  ttf  is  deposited  over  a 
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Fig.  1.  Trilayered  dielectric  structure  with  sources  in  the  cover. 


substrate  region  (><—()  which  is  diaracteiized  by  index  of 
refraction  n,.  The  region  (^  >  0)  is  the  cover  with  refractive 
index  in  which  electric  current  density  J,  maintaining  electro¬ 
magnetic  fields  in  all  three  regions,  is  immersed.  All  media  are 
understood  to  possess  limitingly  small  dissipation.  Although  the 
ensuing  analysis  may  be  generalized  for  a  structure  having  any 
number  of  dielectric  layers  with  embedded  currents,  the  situation 
in  Fig.  1  serves  for  the  purpose  of  illustration,  and  provides  a 
useful  model  for  the  background  of  practical  electronic  and 
optical  integrated  circuits. 

Subsequent  analysis  assumes:  i)  time  harmonic  («-'*")  depen¬ 
dence  of  the  solutions  to  Maxwell’s  equations  and  ii)  all  integrals 
with  unspecified  limits  span  the  entire  space.  The  Hertzian  poten¬ 
tial  subject  to  the  Lorentz  gauge  satisfies  the  Helmholtz  equation 

(  V ^  +  A:?)  n,  -  -  J/juf,  ( 1) 

in  each  region  (i“s,f,c  for  substrate,  film,  cover).  Formal 
operation  on  (1)  with  the  2-D  Fourier  transform 

F{]-jf{}e-'^'dxd2  (2) 

where  \  -  x(  +  rf ,  reduces  (1)  to  the  ordinary  differential  equa¬ 
tion 

y)  ~  y)/ju€^  (3) 

where  v  — F(n},  y  — F{7},  and  pf  —  -  kf.  Solution  of 

(3)  is  elementary,  and  may  be  written  as  a  sum  of  primary 
scattered  parts.  Thus  this  decomposition  is 

v(\.  y)  -  sJ  J  y,r')-^  diA 
\Jy  ju€,  I 

+  +  (4) 

where  g'’(\;  y,r')  -  V2pr,  and  i*  *  Kronecker 

delta.  The  coefficients  fV,^  are  determined  by  satisfying  the 
appropriate  boundary  conditions  [1]  across  the  dielectric  inter¬ 
faces  and  as  y-*  ±oo. 
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Inversion  of  the  tnnsfonn-dtHnain  potentials  yields  the  solu¬ 
tion  to  (1)  with  the  potential  in  the  cooer  region  given  by 

r-  A'-') 

-h/<?'(r|r')~dr.  (5) 

Jy  jci(c 

The  reflected  dyad  G'{r\r')  has  scalar  components  G^(r\r') 
represented  by  2-D  Sominerfeld  integrals  of  the  generic  form 

-^//c^(X)eV^  (6) 

Each  of  the  coefficients  C^(X)  is  a  well-behaved  function  of  X 
in  the  entire  plane;  hence  derivatives  of  G'  may  be  obtained 
by  formally  differentiating  under  the  qiectral  integral.  Special 
attention  is  required  in  determining  derivatives  of  the  principal 
part  of  n. 

The  spectral  integral  mi  the  right  side  of  (S)  represents  the 
primary  part  11'  of  the  Hertzian  potential  It  is  shown  in  the 
Appendix  that  under  the  assunqition  that  /  and  V-J  are  con¬ 
tinuous  and  of  compact  support  in  V,  derivatives  up  to  second 
order  of  II'  may  be  obtained  by  formally  differentiating  under 
the  spectral  inte^raL  Therefore, 

'?vn'(0-^//vy 

le'^'lf  g'(X;;-.r')4^dr 

V  Vy'<y  \) 

where  the  spatial  integration  has  been  split  into  regions  in  which 
g'  is  continuously  differentiable.  Tangential  derivatives  (i.e., 
derivatives  with  respect  to  x  and  z)  of  the  bracketed  term  in  (7b) 
operate  only  on  e^^  '.  However,  performing  the  derivatives  with 
respect  to  y  demands  additional  considerations.  Appropriate  use 
of  Leibnitz’s  rule  (4,  pp.  321-325]  for  differentiation  under  the 
integral  sign  reveals  that 

3  (  [  A'O 

7-  /  g'(X;^,r')-J^dK' 

9yVy‘<y 

r  A'')  \ 

+  / 

Jy->y  JUt^  j 

r  d  /(O 

“/  7- - dV 

Jy<ydy 

{  d  J{r') 

+  /  -^g^(^\y,r')—^dV' 

Jy>ydy 


f-y.iy'-y) 

lA'O 

^Pc 

1  V 

3 

1 

/•  3  J(r') 

^y'<ydy 

t  3  J(r') 

Jy’>y3y  y«<, 

A  subsequent  differentiation  with  reflect  to  y  yields 

3^  (  t  A'O 

XT  /  «'(X;y.r')-^dr 

3y  Vy  <y  J<^*c 

t  A'')  \ 

+/  g^{\;y,r’)—dy\ 

Jy‘>y  jut,  ) 

t  3^  J{r') 

-/  -^g'0-\yy)-^dV' 

^y<y3y  Jf^fc 

t  3^  J(r') 

+  /  7-2«'(X:>,r')-^dr 

Jy’>ydy^  jut. 


l\- 

g-y^y-y') 

3 

g-pAy'-y) 

^Pc 

3y 

A'') 

dx’dz’ 

Me  j  . 

t  3^  ,  /(r') 

Jy<y3y^  jut, 

[  3^  ,  J{r’) 

+/  T^s'0^->yy)—^dv’ 

V>rV  y««c 

-It  — - dx'dz' 

}J  jut. 


so  that  (7b)  becomes 

VVn'(r) 


■//f/f(X;r,r')-^drL^\ 

I'y  jut, 

d^^-LJ{r)/jut,  (8) 


wheie  L  ^  yy  and  the  dyad  g  is  given  by  the  expression 

\  VV  [e^^  ‘'■"'■■’e-'''’'  .  y'>  y- 


The  term  L  J  was  extracted  from  exploitation  of  the  Fourier 
inversion  theorem  [5,  p.  315],  and  is  found  to  correspond  exactly 
with  that  exposed  in  (3]  for  a  “pillbox”  principal  volume.  The 
form  of  g  suggests  that  the  “slice”  exclusion  in  Fig.  2  might  be  a 
more  natural  principal  volume  pertaining  to  L.  This  assertion  is 
verified  in  part  B  of  Section  III. 
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T' 


Fig.  2.  A  "slice''  prindpsl  vdiime  excluding  (be  singulafity  point  of  the 
electric  dyadic  Green's  function;  closed  surface  is  the  boundary  of  the 
slice  volume. 


"slice”  volume  shown  in  Fig.  2.  Starting  with  the  common 
representation  for  the  free-space  Green’s  function  f  (r|r')  - 
—  rX  it  may  be  shown  that  for  a  “slice"  principal 
volume,  the  correction  term  £'(r)  for  the  electric  dyadic  Green’s 
function  for  fleld  points  in  the  source  region  is 

1  , 

£^(r) - lun  /  V’H>-\r’)h'  J{r')dS'  (13) 

ytix^  «  -  0  •'s, 

where  is  shown  in  Fig.  2.  The  correction  term  above  is  now 
shown  to  be  equivalent  to  the  correction  term  corresponding  to  a 
pillbox  principal  volume.  The  surface  integral  term  in  (13)  b  split 
into  integrations  over  5,  and  (planes  at  y±S,  respectively). 
This  yields  for  (13) 

E^ir) —  hm  l~f  V'Hr\r%-y*a 
yu<<.a-o(  Jsy 

■J^ix\y  +  S,z')dx'dz' 


III.  Electmc  Dyadic  Green’s  Function 
A.  Development  of  the  Principal  Dyad 

The  electric  Geld  E  is  related  to  the  Hertzian  potential  by 
E  “■  (ikj  +  W-)n.  Using  (8),  the  principal  part  of  the  field  may 
be  written  as 

£'(  r)  -  -  f  ( 4i'( X:  r.  r')  •/(r')  drj  d^\ 

(10) 

where  the  dyad  g’  '“g/k^  +(l/4ir*)/g'’e^^  '. 

Equation  (10)  is  a  useM  expression  for  the  principal  part  of 
the  electric  field  due  to  the  simple  nature  of  the  integrand 
appearing  in  the  volume  integral  However,  (10)  is  not  written  in 
the  standard  form  as  a  volume  integration  of  Ihe  scalar  product 
of  a  Greeks  dyad  with  the  electric  current  density.  The  depolariz¬ 
ing  dyad  L  has  manifested  itself  naturally.  Recognizing  that  the 
corresponding  "principal  volume”  is  a  pillbox  [3]  yields  the 
standard  form  for  the  electric  field: 


£'(r)  -  -  jupo  lim  /  G'(r,  r')  /(r')  dV'-L  ■J(r)/ju€^ 

v  —  O'^y-v 

(11) 

where  u  is  a  pillbox  excluding  the  singularity  of  G*  at  r  and  C' 
is  given  by 


G'(r,r') 


(/-t- vv//c^) 


g-pAy-y') 


2(2ir)^/>, 

g-pAy'-y) 

2(2ir)^/>c 


d^X, 

d^X, 


y'<y 

y'>y. 

(12) 


B.  Equivalence  of  Principal  Volumes 

The  principal  volume  t;  in  (11)  was  identified  to  be  a  pillbox  as 
tabulated  in  [3].  A  more  useful,  and  equivalent,  exclusion  is  the 


As  S  -►  0,  Sx  approaches  and  Jy{x',  y  +  S,z')  approaches 
Jy(x‘,  y  —  6,z')  due  to  the  smoothness  of  the  boundary  of  V  and 
the  continuity  of  /  at  y'  •  y.  Thus,  (14)  simplifies  to 

£^(r)-—  lim  /  vV(rlr')|;::p*/,(x',y,r')dx'dr'  (15) 
Jb)t^  S  -  0  ■'s  '  ^ 

where  S  extends  over  the  x’-z'  plane.  Expressing  vV  in  Carte¬ 
sian  form  as 

V>(rlr')-(-yA:,-l/£)^^ 

•(x(x'-  x)  +  y{y’-y)  +  f(r'-  r)] 
where  £  —  |r  —  r'j,  it  is  found  that 

V'+(rlr')i;;:>:*  -  (  -  jk,-l/R)^^y2S  (16) 

where  £,  —  ((x  —  x*)^ -t-(r  —  r')^ Substitution  of  (16) 
into  (15)  yields 

£'(r)  -  —  lim  (  y/2«(  -  jk,  -  l/R,) 

JU€^  S  0  (  ■'S 

The  integral  in  (17)  may  be  decomposed  into  the  sum  of 
integrals  over  5 -  Q  and  c;.  Q  is  a  circle  centered  at  (x,  z)  with 
radius  p.  As  8  -*  0,  integration  over  S  ~  C,  vanishes.  If  v  is 
chosen  sufficiently  small,  then  /,(x',  y,  z')  *  /^(r)  and  = 

1  so  that  (17)  becomes 

£<-(^)___ - j>/(r)  Urn  (sf  — — dx'dz]  (18a) 

ju(^  '  t^o\  Jc.  2nR]  f  ''  ' 


- - K('’)  Jim 

juc^  «  -  0 


1  if*  J^c 

hm|fi/  - - - pdp 
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when  Tf  -  (p?  +  In  going  from  (18a)  to  (18b),  integration 
over  C,  has  been  transformed  to  polar  coordinates.  Performing 
the  angular  integratkm  is  trivial  and  yields  (18c).  Noting  that  the 
integrand  in  (18c)  is  a  perfect  differential,  the  term  in  braces 
becomes 

-«{  yikJln(i-^  +  8*)‘'^-ln«] 


where  f{v-/}  is  the  Fourier  transform  of  V  /  as  deflned  by 

(2). 

Next,  use  p,  “  p,  +  jpi,  where  p,  and  p,  are  the  real  and 
imaginary  parts  of  p„  respectively.  The  exponential  is 

of  constant  sign  for  all  y.  By  the  generalized  first  mean  value 
theorem  for  integrals  [4,  p.  117],  the  right  side  of  (21)  may  be 
written  as 


+  [l/«-(r*  + 
-1  (as«-0). 


d^\ 


Finally,  substitution  of  (19)  into  (18c)  yields 

y«<c 

which  is  precisely  the  same  correction  term  appearing  in  (11). 
Therefore  integration  excluding  a  “slice”  principal  volume  is 
equivalent  to  a  pillbox  exclusion,  and  evaluation  of  (11)  may  use 
either  of  these  volumes. 


L 

where  p„u„  <  (n.^)  <  (-^“0  for 

The  spatiid  integration  in  (22)  is  trivial  and  leads  to 


rv.  Conclusions 

In  the  study  of  layered-media  electromagnetics,  Sommerfeld 
integrals  are  used  to  represent  scalar  components  of  the  electric 
dyadic  Green’s  function.  The  principal  part  of  the  electric  Held 
may  be  written  as  a  sum  of  a  spectral  integral  along  with  a 
correction  term  that  appears  naturally.  The  spectral  integration 
may  be  replaced  with  the  more  standard  volume  integral  as  in 
(11).  Recognition  of  the  depolarizing  dyad,  which  has  manifested 
itself  innately,  identifies  the  appropriate  principal  volume  (a 
pillbox).  An  equivalent  excluding  region  to  the  pillbox  is  sug¬ 
gested  to  be  useful  in  practice  due  to  its  simple  form. 

Appendix 

It  is  now  shown  that  the  differentiation  under  the  specual 
integral  as  in  (7)  is  a  legitimate  operation.  Without  loss  of 
generality,  justifying  this  interchange  of  operations  for  the  follow¬ 
ing  is  sufHcient: 

//  W  e^^  'f  g»i\-,y.r')J{r')dy  d^\.  (Al) 

•' •'k  >  k,  I  •' y 

In  (AI),  k,  is  the  real  part  of  Evaluation  of  is  made  on  the 
Riemann  sheet  with  Re  ( p, )  >  0. 

Assuming  that  J  and  V-/  are  continuous  and  have  compact 
support  in  V,  use  of  the  vector  identity  V-  (9^4 )  —  99  /4  -t-  V9  d 
along  with  the  divergence  theorem  on  (Al)  yields 

V'  Jir’)g''i\;y,r’)dy'  d^X 

■'■'k>k,  [  ■  ■'y 

"If  ^  e'^^'lFiv'-Jir'))^ - dy'  d^X  (A2) 

•'  •'X  >  *,  •'  iPc 


/jf>^V  e^’’^Re(F{v  /(x'.,.r')} 


where  g>{pr',y)“  (2- Since  v  J 's 
continuous  and  of  compact  support  in  K,  v  J^L}  (i.e.,  the 
space  of  square  integrable  functions).  In  particular,  for  each  y, 
V-/  is  an  L?  function  in  the  i-z  plane.  Using  a  standard 
theorem  from  Fourier  transform  theory  [S,  pp.  310-313],  the  2-D 
Fourier  transform  of  V-/  is  an  L}  function  in  the  plane. 
Thus,  F[v-J)  “0(X"*”*)  as  (X-»oo,  oO).  The  integrand  in 
(A4)  is  dominated  in  magnitude  by  a  function  which  is  indepen¬ 
dent  of  r  and  0(X”^"*).  The  Weierstrauss  M  test  [4,  p.  470] 
guarantees  that  the  integral  in  (A4)  converges  uniformly,  whereby 
a  standard  theorem  from  advanced  calculus  [4,  p.  474]  justifies 
the  interchange  of  differentiation  and  spectral  integration. 
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In  the  study  of  layered  environments  in  electronic  and 
optical  integrated  circuits,  integral-operator  formulations  are 
increasingly  exploited.  Construction  of  such  integral  operators 
requires  an  appropriate  Green's  function.  In  this  paper,  we  pre¬ 
sent  the  development  of  a  general  electric  dyadic  Green's  function 
for  layered  structures  with  magnetic  as  well  as  dielectric  con¬ 
trast.  The  general  properties  of  that  Green's  dyad  are  studied. 

In  the  tri-layered  environment,  a  film  of  thickness  t  with 
constitutive  parameters  is  deposited  over  a  substrate 

(y<-t)  described  by  The  region  {y>0)  is  the  cover  with 

parameters  (c  ,u  ).  The  arbitrarily  directed  electric  current  ^ 
c  c 

is  immersed  in  the  cover  or  film  region.  By  Sommerfeld's  classi¬ 
cal  method,  a  Hertzian-potential  dyadic  Green's  function  is  found 
whih  satisfies  general  boundary  conditions.  The  electric  Green's 
dyad  is  expressed  in  terms  of  its  Hertzian  counterpart  as 

^(rlr')  =  (k^  +  W.)g(rlr')  +  C6(r-r') 
c 

?(rl?')  =  fGP(rlr')  +  ^{v,P) 

where  ?  is  the  Hertzian-potential  Green's  dyad,  L  is  a  depolariz- 

ing  dyad,  G*^  is  the  principal  Green's  function,  and  ?  is  a  re¬ 
flected  Green's  dyad.  The  scalar  components  of  the  reflected  dyad 
are  obtained  in  2-D  spectral  integral  representations.  The  dyadic 
Green's  functions  are  reducible  in  certain  special  cases,  such  as 
that  of  a  conducting  substrate  layer. 

The  source-point  singularity  of  the  electric  Green's  dyad  is 

studied,  and  depolarizing  dyad  L  is  identified.  It  is  demonstrat¬ 
ed  that  the  electric  Green's  dyad  for  layered  media  of  general 
contrast  is  reciprocal.  Pole  singularities  of  implicated  reflect¬ 
ion  coefficients  lead  to  general  discrete  surface-wave-mode  eigen¬ 
value  equations. 

The  electric  dyadic  Green's  function  obtained  is  applied  to 
the  microstrip  transmission-line  structure  with  magnetic  (ferrite) 
film  and  conductor  substrate.  The  EFIE's  based  upon  that  Green's 
function  are  solved  by  the  HoH  and  by  a  perturbation  approxima¬ 
tion.  Numerical  results  are  obtained  and  compared  with  those  of  a 
structure  having  a  purely  dielectric  surround. 
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DYADIC  GREEN’S  FUNCTIONS  FOR  EM  FIELDS 
IN  THE  LAYERED  PC/IC  ENVIRONMENT  -- 
REPRESENTATION,  PROPERTIES  AND  APPLICATIONS 

Dennis  P.  Nyquist 

Department  of  Electrical  Engineering 
Michigan  State  University 
East  Lansing,  Michigan  48824 


The  study  of  EM  interactions  among  electronic/optical  devices  in  the 
layered  PC/IC  environment  has  become  a  contemporary  topic  in  electromagnetics 
research.  Interest  is  in  microstrip  or  dielectric-waveguide  based  circuits 
operating  at  microwave,  millimeter,  or  optical  wavelengths.  Such  circuits  are 
typically  located  adjacent  to  a  planar,  tri-layered  substrate/film/cover 
surround.  In  microwave  and  most  mm-wave  circuits,  the  dielectric  substrate 
layer  is  replaced  by  a  good  conductor.  This  paper  exposes  dyadic  Green's 
functions  useful  or  the  integral-operator  description  of  such  circuits,  as 
well  as  certain  oi  their  general  properties  and  applications. 

Dyadic  Green's  functions  for  the  Hertz  potential  and  electric  field 
maintained  by  electric  currents  immersed  in  a  tri-layered  substrate/film/cover 
surround  are  obtained  in  Sommerfeld  integral  representation;  they  are  cast 
into  a  form  useful  for  the  description  of  the  EM  field  in  the  PC/IC  environ¬ 
ment.  Arbitrary  contrast  of  dielectric/magnetic  properties  among  the  several 
linear,  isotropic,  homogeneous  media  in  the  tri-layered  surround  is  accomo¬ 
dated.  A  substrate  with  negative-infinite  imaginary  permittivity  provides  the 
special  case  of  a  conducting  layer.  Alternative  representations,  arising  from 
real-line  Fourier  inversion  and  its  subsequent  conversion  to  a  singularity 
expansion  involving  pole  and  branch-point  singularities,  are  exposed. 

Several  general  properties  of  the  Green's  dyads  are  identified.  Pole 
singularities  of  integrands  in  the  dyad  representations  correspond  to  TE/TM 
surface  waves  supported  by  the  asymmetric  substrate/film/cover  surround  with 
appropriate  dielectric/magnetic  contrast  among  layers.  The  electric  Green's 
dyad  is  demonstrated  to  be  reciprocal,  and  its  source-point  singularity  is 
exposed.  The  latter  singularity  leads  to  an  innate  contribution  to  the 
electric  field  at  a  source  point  which  is  absent  in  conventional  eigenfunction 
expansions  of  that  field.  A  mixed-potential  formulation  identifies  a  scalar 
potential  maintained  by  charge  immersed  in  the  tri-layered  environment. 

The  Green's  dyads  are  applied  to  obtain  an  integral-operator  description 
for  electronic/optical  components  in  PC/IC  configurations.  Application  to 
open  waveguides  leads  to  quantification  of  the  complete  propagation-mode 
(including  radiation)  spectrum  as  well  as  an  excitation  theory  for  discrete 
nodes.  A  coupled-mode  perturbation  theory  for  systems  of  adjacent  open  wave¬ 
guides  is  exposed.  The  integral-operator  formulation  also  provides  a  descrip¬ 
tion  for  coupling  of  patch  antennas  and  other  devices  among  themselves  and 
associated  transmission  systems.  Typical  numerical  results  will  be  presented. 
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A  NOTE  ON  THE  MIXED  POTENTIAL  REPRESENTATION  OF 
ELECTRIC  FIELDS  IN  LAYERED  MEDIA 

M.J.  Cloud  and  D.P.  Nyquist 
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ABSTRACT 


A  mixed  potential  formulation  is  given  for  electric  fields  in 
layered  environments.  Contributions  to  the  field  from  charges  are 
identified  explicitly  through  a  scalar  Green's  function  for  layered  media. 
The  outcome  is  a  computationally  expedient  Sommerfeld  integral 
representation . 
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I.  Introduction 

The  study  of  electric  fields  due  to  surface  currents  in  millimeter- 
wave  integrated  circuits  [1,2]  brings  to  light  some  facts  about  the 
alternative  representation  of  Hertz  potentials.  Despite  their  apparent 
simplicity,  these  observations  have  not  appeared  previously  in  this  form. 

The  points  discussed  in  this  paper  bear  directly  upon  the  divergent  spectral 
integrations  which  are  offered  on  several  occasions  [3,4]  in  the  recent 
literature;  it  is  hoped  that  ultimately  they  will  find  application  in 
avoiding  comparatively  awkward  formulations. 

II.  Hertz  Potential  Green's  Dyad  Formulation 

Consider  the  configuration  of  layered  dielectric  media  over  a 
conducting  half-space  as  shown  in  Figure  1.  The  electric  field  E(r)  in  the 
cover,  maintained  by  surface  currents  embedded  in  that  same  layer, 
decomposes  linearly  into  two  parts  as  E(r)  -  E^(r)  +  E’'’(r) .  The  fields  of 
the  right  member  may  be  termed  the  primary  and  reflected  components.  A 
Hertz  potential  representation  of  E  based  upon  this  decomposition  is  given 
by  Bagby  and  Nyquist  [1,2]  as  follows; 

l(r)  -  -t-  W.)  |gG(?,?')-[K(?')/jwc]dS'  (1) 

where  K(r)  describes  source  currents  on  surface  S,  and  G  -  IG^  +  G  is  the 
decomposition  of  the  Hertz  potential  Green's  dyad  into  primary  and  reflected 
components.  The  dyad  scalar  components  are  given  as  double  spectral 
(Sommerfeld-type)  integrals: 

o  -  -  rr  exp[jA.(r  -  r')]exp[-p  |y  -  y' 1 ]  „ 

GP(r.r')  -  - - - 2 -  d^X  (2) 

ii  2i2n)\ 


23 


G^(?.r') 


R  (A) 
R„(A) 
C  (A)J 


exp[jA«(r  -  r*)]exp[-p^|y  -  y’ | ] 

2(2x)V 


d^A 


(3) 


AAA  A  A  AAA 

where  -  xG^  x  +  y[(3G^/ax)  x  +  G^y  +  (aG^/3z)2]  +  z  G^  z.  The  R^,  R^, 

and  C  are  reflection  and  coupling  coefficients  detailed  in  [1,2].  Note  that 

2  2  1/2  '' 

p^  -  [A  -  k^]  '  is  a  wavenumber  parameter,  A  -  xA^  +  zA^  is  a  vector 

2 

spatial  frequency,  and  d  A  -  dA^dA^. 

The  first  term  (with  factor)  in  E(r)  is  proportional  to  the 
vector  Hertz  potential  contribution  to  E  (produced  directly  by  currents  K) , 
while  the  second  term  is  subsequently  treated  as  the  negative  gradient  of 
scalar  potential.  In  the  present  formulation,  and  in  many  similar  ones 
[3,4],  it  is  found  that  boldly  exchanging  the  order  of  differentiation 
(i.e.,  the  TV*  operation)  with  integration  too  many  times  leads  to 
convergence  problems  with  the  resulting  inner-nested  spectral  integrals. 

Care  will  be  taken  in  the  present  development  to  allow  only  one  such 
interchange;  moreover,  the  impact  of  the  derivative  on  the  Green's  function 
is  subsequently  minimized  by  a  technique  of  integration-by-parts. 

Figure  2  illustrates  the  geometry  of  surface  S,  supporting  currents 
of  integrated  electronics.  Assume  that  K(r)  is  a  continuous  function 
everywhere  on  S,  "xcept  possibly  at  the  boundary  of  subregion  S^.  The 
exception  is  made  because,  although  one  physically  expects  continuous 
currents,  discontinuities  could  be  introduced  mathematically  by  subsectional 
basis  functions  in  numerical  solutions.  and  are  auxiliary  boundary 

contours  just  outside  and  inside  of  C  .  Take  Cq  to  be  the  outer  contour  of 

A 

S.  Let  S'  -  S*Spi  with  n'  the  outward  normal  unit  vector  to  S' .  It  is 
important  to  note  here  that  principal  Green's  function  may  in  fact  be 
represented  either  in  Sommerfeld  integral  form,  or  simply  as  exp(  -  jkR)/4;rR; 
the  latter  form  points  out  explicitly  the  nature  of  its  source-point 
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(R  -  |r  -  r' I  -  0)  singularity.  This  singularity,  unless  handled  properly, 
invalidates  a  subsequent  application  of  the  divergence  theorem.  To  maintain 
complete  rigor  a  small  area  should  be  excluded  from  S,  at  the  location  r,  to 
preclude  the  coincidence  of  r  and  r'  during  the  spatial  integration  in  (1) . 
However,  at  points  where  K  is  continuous  there  is  no  contribution  from 
either  the  excluded  area  or  its  boundary  contour.  Therefore,  the  exclusion 
process  is  not  detailed  here.  Points  of  surface  current  discontinuity  are 
handled  in  the  following  section. 

Ill .  Primary  Component  of  Electric  Field 

Consider  first  the  primary  component  of  E.  The  electric  scalar 
potential  contribution  to  E^  is  where 

-  -<l/J«OV-JgGPKdS' 

-  -(l/jwO[Js.V*(GPK)dS'  +  L  V.(GpK)dS'].  (4) 

J  J  p 

This  expression  is  subsequently  cast  into  a  form  which  explicitly  exposes 
its  electric  charge  sources.  Although  the  method  is  conventional,  it  is 
included  here  to  guide  the  identification  of  similar  contributions  to  the 
reflected  potential.  The  procedure  is  less  obvious  in  the  latter  case. 

Exchanging  the  divergence  and  spatial  integration  prompts  vector 
manipulation  of  the  integrand  as  follows: 

V*(GPit)  -  G^V*!^  +  VgP.K  -  -V'*(GPi^)  +  G^V' -i!.  (5) 
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Into  the  last  term  a  substitution  from  the  surface  continuity  equation, 
7<K  -  is  made.  The  scalar  potential  becomes 


♦P  -  (l/j«c)[L,7'*(GP2)dS'  +  fg  7'.(GPK)dS'] 
J  J  p 


+  (lA)Jg(GPa)dS' 


(6) 


to  which  the  2-D  divergence  theorem  is  applied  to  give 


-  (l/0[(l/jw)  |^^GP(n'.it)dl' 

+  (1/jw)  GP(n'-(^  -  K’)dl' 

+  JgGPadS'] 

-  (l/0[  cPpj^pdl'  +  JgGPadS']  (7) 


where  K  and  K  are  values  of  surface  current  on  and  ,  respectively. 
Thus  *P  is  written  in  terms  of  line  and  surface  charges,  while  the  original 
derivatives  have  been  effectively  removed  from  gP. 


IV.  Reflected  Component 

For  the  reflected  potential. 


* 


r 


(-l/j«€)7. 


K  dS' 


(8) 
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one  may  follow  the  same  general  sequence  of  steps  as  above,  but  this  time 

exploiting  the  specific  form  of  G  .  Following  the  procedure  beginning  with 

equation  (5),  and  defining  a  scalar  Green's  function  g’"  -  G^  +  dG^/dy, 

c  o 

v.(g’^.k)  -  v.{xG^K^  -  J[(aG^/ax)K^  +  (aG^/a2)K^]  +  ^g^k^) 

-  )aGj/3x)K^  +  (aG^/az)K^ 


+  (a/ay)f(aG^/ax)K^  +  (aG^/3z)K^] 

-  ^  [\^(8Gl/8y)]-K 

-  (VG^)*K  -  (.7'G^)«K 


(9) 


Carrying  the  analogous  steps  to  completion,  one  finally  arrives  at 

-  (l/£)[  I^^G^'p^Qdl'  +  G^'pj^pdl'  +  JgG^'adS']  (10) 

Therefore  the  reflected  scalar  potential  is  expressible  in  terms  of 
effective  charges  weighted  by  an  appropriate  Green’s  function.  It  is  noted 
that  convenient  form  of  the  scalar  reflected  Green's  function  arose  from  the 
specific  form  of  the  Hertz  potential  Green's  dyad  cited  earlier.  Since  the 
conversion  process  has  effectively  removed  derivatives  from  this  function, 
the  new  formulation  is  expected  to  be  computationally  efficient.  Finally, 
the  ^  field  in  the  cover  region  is  exp  ressed  as 

JgGodS']  (11) 


i  -  V(l/£)(|^^Gp^Qdl'  +  Gpj^pdl' 
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where  G  -  -i-  is  a  Sonmerfeld  integral  representation  of  the  scalar 

potential  Green's  function  for  layered  media.  The  gradient  operator  may  (at 
interior  points)  be  exchanged  with  the  spectral  Integrals  without  rendering 
those  Sommerfeld  integrals  non- convergent. 

V.  Summary  and  Conclusions 

The  preceding  sections  illustrate  a  procedure  for  avoiding  the 
needless  imposition  of  derivatives  onto  the  Green's  function,  thus  avoiding 
convergence  problems.  In  the  process  of  converting  to  a  less  singular 
formulation,  unknown  surface  and  line  charges  are  introduced  explicitly  into 
the  problem.  Presumably,  for  many  applied  problems,  these  charge  functions 
could  be  expanded  in  suitable  moment  method  basis  sets  along  with  the 
original  surface  currents. 
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CONDUCTOR 


Figure  1.  Configuration  of  Layered  Media 
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4.2  Mathematical  methods — analytical  and  numerical 


Efforts  on  this  topic  were  devoted  to  the  formulation  and  computation  of 
Sommerfeld-integral  type  representations  for  relevant  Green’s  function  ker¬ 
nels.  Analytical/numerical  methods  to  accomodate  the  slowly-convergent  (and 
sometimes  apparently  divergent)  behavior  of  those  representations,  arising 
when  the  source-point  singularity  is  implicated,  were  developed.  These  init¬ 
ial  preparatory  investigations  were  necessary  to  acquire  confidence  in  the 
construction,  computation,  and  use  of  the  dyadic  kernels  described  in  Section 
4.1.  Studies  on  the  computation  of  Green’s  function  kernels  are  reported  in 
[1,8,12],  while  properties  of  the  source-point  singularity  are  addressed  in 
[8]. 
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ACCURATE  EVALUATION  OF  SOMMERFELD  INTEGRALS 

USING  THE  FAST  FOURIER  TRANSFORM 

B.C.  Drachman 

Department  of  Mathematics 

D.P.  Nyquist  and  M.J.  Cloud 

Department  of  Electrical  Engineering  and 

Systems  Science 

Michigan  State  University 

East  Lansing,  MI  48824 

Application  of  the  conventional  Hertzian-poten- 
tial-based  formulation  of  layered  media  electromag¬ 
netics  necessitates  evaluation  of  Sommerf eld-type 
integrals.  An  efficient  and  accurate  method  of 
computing  these  integrals  is  presented. 

The  Fast  Fourier  Transform  (FFT)  may  be  used  in 
a  straight-forward  fashion  to  obtain  accurate  inte¬ 
grations  of  periodic  functions.  It  is  explained 
how,  in  these  cases,  the  numerical  integration  meth¬ 
od  corresponds  to  the  Corrected  Trapezoidal  rule 
(with  error  proportional  to  the  fourth  power  of  the 
step  size).  Attempts  to  use  the  FFT  in  such  a 
fashion  on  aperiodic  functions  result  in  significant 
errors  because,  in  these  cases,  the  integration 
method  corresponds  to  the  Rectangular  rule  (error 
proportional  to  the  first  power  of  the  step  size). 
Thus  it  is  often  believed  that  the  FFT  is  a  quick 
but  inaccurate  way  to  compute  an  integral.  When 
used  in  conjunction  with  Simpson's  rule,  however, 
the  FFT  gives  extremely  accurate  results  for  aperi¬ 
odic  functions.  An  error  analysis  suggests  a  tech¬ 
nique  for  integral  truncation  and  choice  of  parti¬ 
tion  density.  Examples  are  given  in  which  this 
method  is  applied  to  the  numerical  evaluation  of 
Sommerf eld  integrals. 

The  method  exploits  the  efficiency  of  the  FFT 
algorithm,  while  deriving  additional  benefit  from 
the  fact  that  a  single  call  of  the  FFT  algorithm 
returns  sample  values  of  the  Sommerfeld  integral  at 
many  spatial  locations.  Moreover,  this  technique 
should  be  well  adapted  for  use  with  the  point¬ 
matching  Galerkin's  method  of  moments  which  requires 
values  of  the  Sommerfeld  integrals  at  a  set  of 
locations  on  a  structure. 
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ON  SOMMERFELD-INTEGRAL  ELECTRIC  FIELD  KERNELS  FOR 
HICROSTRIP-BASED  CIRCUITS 

D.P.  Nyquist,  M.S.  Viola,  M.J.  Cloud  and  M.  Havrilla 
Department  of  Electrical  Engineering 
Michigan  State  University 
East  Lansing,  MI  48824 

Microstrip  circuits  are  fabricated  over  the  film/cover 
interface  of  a  planar,  layered  conductor/film/cover  surround 
environment.  In  an  electric-field  integral  equation  descrip¬ 
tion  of  such  circuits,  the  electric  field  tangential  to  metal¬ 
lic  microstrip  conductors,  maintained  by  surface  currents  on 
those  conductors,  is  required.  A  spectral  approach  through 
Hertz  potentials  leads  to  sommerf eld- type  integral  representa¬ 
tions  for  required  Green's  function  kernels. 

Let  the  microstrip  conductors  reside  at  y»0^,  adjacent  to 
the  z-z  film/ cover  interface  plane.  The  field  e(x,c)  at  the 
conductor  surface,  maintained  by  surface  current  k(x,^)  on 
that  conductor  (both  in  the  axial  Fourier-transform  domain; 
transform  variable  (;),  is  then 

e(x,?)  =  (k^+VV*)  J  ^  g^(x,ylx',0)*k(x',i;)dx' 

where  k  =n  k^  (n  =cover  index) ,  V  =  and  g.  is  a  Hertz 

potential  Green's  dyad  in  Sommerf eld^ntegral  representation. 
It  is  convenient  to  pass  the  spatial  integral  and  the  differ¬ 
ential  operator,  in  turn,  through  the  spectral  integral  over 
since  the  former  operations  can  all  be  performed  in  closed 
form.  This  procedure  is  valid  provided  the  details  |source- 
point  singularity)  are  handled  carefully  and  the  y-»-0  limit  is 
performed  lastly.  In  much  of  the  recent  literature,  the  lat¬ 
ter  limit  is  exchanged  with  the  spectral  integral;  this  inval¬ 
idates  the  differential-operator  exchange  and  leads  to  a  di¬ 
vergent  spectral  integration.  Consequently,  uncertainty  has 
arisen  in  the  choice  of  appropriate  expansion  and  testing 
functions  for  use  in  MoM  solutions  to  associated  EFIE's. 

In  fact,  pulse  basis  functions  and  foint  matching  lead  to 
convergent  spectral  integrals  if  the  y-^  limit  is  performed 
appropriately.  MoM  matrix  elements  are  represented  by 


-p^y 


where  V  .  is  a  complicated  function20f2^. (carying  surface-wave 
poles  or  the  surround)  and  p  =(C  .  The  most  slowly 

convergent  case  is  o=8=x,  whire  for  large  arguments, 

leading  to  an  undamped  oscillator^^ntegrand  when  y^O.  It  is 
demonstrated,  however,  that  the  y-»-0  limit  does  indeed  exist, 
and  the  convergence  rates  of  the  various  I  ^  are  compared. 
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Accurate  Evaluation  of  Sommerfeld  Integrals  Using  the 
Fast  Fourier  Transform 

B.  DRACHMAN*t**,  M.  CLOUD*,  AND  D.  NYQUIST* 


Abstract  -  Accurate  numerical  integration  of  Sommerfeld-type  integrals  by 
the  Fast  Fourier  Transform  algorithm  is  considered.  Several  examples 
illustrate  in  the  process. 


I.  INTRODUCTION 

The  purpose  of  this  note  is  to  observe  that  the  FFT  combines  naturally 
with  Simpson’s  rule  for  Sommerfeld-type  integral  computation.  Standard 
references,  such  as  {1,2,3],  discuss  the  FFT  for  periodic  functions  only; 
however,  even  relatively  simple  Sommerfeld  integrals  have  aperiodic 
integrands. 

A  principal  advantage  of  using  the  FFT  is  that  a  single  subroutine  call 
yields  a  set  of  sample  values  of  an  integral  (i.e.,  the  integral  for  various 
values  of  an  integrand  parameter).  Such  samples  could  be  useful  in 


*This  work  partially  cupported  by  Office  of  Naval  Research  under 
Contract  N0014-86-K-0609. 

**Partially  supported  by  DARPA. 

B.  Drachman  is  with  the  Department  of  Mathematics,  M.  Cloud  and 
D.  Nyquist  are  with  the  Department  of  Electrical  Engineering,  Michigan 
State  University,  East  Lansing,  Michigan  48824. 


35 


themselves;  for  example,  the  parameter  mi^ht  be  a  spatial  coordinate  variable 
in  a  moment  method  solution.  In  other  applications  Sommerfeld  integrals 
represent  Green’s  functions  nested  within  other  spatial  integrals.  Samples 
from  the  PFT  might  be  useful  in  approximating  the  outer-nested  integral. 


II.  SUMMARY  OF  NEWTON-COTES  NUMERICAL  INTEGRATION  METHODS 
To  estimate  an  integral  such  as 

I  =  J  f(t)dt  (1) 

a 

suppose  [a,b]  is  partitioned  into  N  subintervals  [ti>ti+i]  with 
ti  =  a  +  hi,  i  =  0,1,2,.. .,N-1.  The  step  size  h  is  {b-a)/N.  Table  1 
summarizes  several  well-known  integration  rules,  provided  f(t)  is  sufficiently 
smooth.  Note  that  rj  is  some  point  in  (a,b},  and  N  is  even  for  Simpson’s 
rule.  These  rules  are  summarized  here  for  reference  in  later  sections. 


III.  NUMERICAL  INTEGRATION  WITH  THE  FFT 

The  FFT  algorithm  computes  the  discrete  Fourier  transform 
N-1 

=  I  exp(j2TTkn/N);  k  =  0,1,2 . N-1  (2) 

n=0 


of  the  sequence  A,,...,A[^.  This  may  be  viewed  as  numerical  integration  to 
approximate 

I(x)  =  J  g(t)exp( jxt)dt  (3) 

0 
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as  follows.  First,  choose  T  sufficiently  large  so  that  the  improper  integral 
in  (3)  is  adequately  approximated  by 

I(x)  =  j  g(t)exp(jxt)dt  (4) 

0 

for  X  in  the  range  of  interest.  Make  a  change  of  variables,  transforming 
the  interval  of  integration  to  [0,27r]: 

I(x)  =  (T/2n)  J  g(tT/27i)exp( jxtT/27T)dt  (5) 

0 

Partition  this  new  interval  into  N  equal  segments,  each  of  width  27r/N.  The 
rectangular  rule  (see  Table  1)  then  yields  the  approximation 

N-1 

I(x,  )  =  h  I  g(nh)exp(j2TTnk/N)  ,  (6) 

^  n=0 

where  h  =  T/N  and  x^  =  2TTk/T.  Comparison  of  equations  (1)  and  (6)  shows 
I(xk)  =  hS^+i,  where  An+i  =  g(nh).  Therefore,  the  FFT  is  directly  invoked 
for  efficient  computation  of  I(x)  by  the  rectangular  rule. 

The  great  accuracy  found  by  those  using  the  FFT  for  periodic  functions 
is  easily  explained.  Observe  (see  Table  1)  that  if  g(t)  is  periodic  and 
sufficiently  smooth  on  [0,T],  the  rectangular  and  corrected  trapezoidal  rules 
are  identical.  Thus  fourth-order  accuracy  is  obtained.  For  aperiodic  g(t) 
the  method  does  not  correspond  with  corrected  trapezoidal  rule,  but  to 
rectangular  rule,  and  is  therefore  accurate  only  to  first-order. 

For  g(t)  aperiodic,  better  accuracy  is  obtained  by  weighting  the 
integrand  sample  values  g(nh)  prior  to  the  FFT  call.  Reference  to  Table  1 
shows  that  the  simple  procedure  of  letting  A,  =  g(0)  +  g(hN),  A,  =  4g(h), 
As  =  2g(2h),  A4  =  4g(3h)  ,  ...,  implements  Simpson’s  rule  within  a  factor 
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of  h/3 


To  demonstrate  the  consequences  of  this  improvement,  consider  the 
following  example:  In  order  to  compute  the  transform 

\  1  f“  -t^/2  jxt 

g(x)  =  J  e  '  e**  dt 

y/2TT  — «» 

=  -L  Real  f*  ^jxt 

V27r  0 

which  we  approximate  by 

f  dt 

V2x  0 


(Take  the  real  part.)  We  choose  T  =  4n  and  N  =  2^®.  The  FFT  computes  gCxj^) 
at  1024  sample  points.  To  check  accuracy,  let  xj^  =  1.  Then 


Since  by  (2),  Sj^+j 


27Tk  _  , 
T  ■ 

N-1 


^2TTjik/N 


so  k  =  2. 


i=0  _ 

S3  corresponds  to  gCx^),  that  is  gCx,)  =  ^  Ss- 

-t^/2 

Except  for  a  constant  factor,  e  is  its  own  tremsform.  Hence  g(l) 

should  be  e"*5.  Direct  FFT  with  unweighted  integrand  sample  values  (i.e. , 
rectangular  rule)  produces  an  answer  correct  to  only  one  significant  figure. 
The  Simpson’s  rule  scheme  gives  13  digit  accuracy  using  a  CDC-750  computer  eind 
the  FFT  routine  from  the  IMSL  library. 


For  another  example,  suppose  sample  values  of  Ko(x)  are  needed,  where 
Ko  is  the  modified  Bessel  function  of  order  0,  a  Green’s  function  for  the 
Helmholtz  wave  equation.  For  example,  the  sample  values  might  be  needed  to 
carry  out  a  numerical  integration  involving  Ko(x)  in  another  problem. 
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Although  there  are  better  methods  for  estimating  Ko(x),  we  use 


jxt 


dt 


(7) 


In  order  to  make  the  integrand  decay  more  quickly,  use  integration  by 
parts  four  times  and  obtain 


2 

3  o 


-f 


8t  -24t  +3  jxt 


dt 


(8) 


so 


K 


/■o = ^  r 


8t  -24t  +3  jxt  .. 

- §72“  ® 


(9) 


with  «  replaced  by  T  =  4Tr  and  N  =  1024,  a  single  FFT  call  with  Simpson’s 
rule  returned  values  Ko(l),...,l(o(10)  to  6  place  accuracy. 


IV.  ERROR  ANALYSIS 

Two  types  of  error  are  introduced  in  the  previous  computations: 
truncation  and  resolution  errors.  First,  consider  the  truncation  error.  The 
error  introduced  by  truncating  an  improper  integral  is  bounded  by  proper 
choice  of  the  stopping  point  T.  Assuming  g(t)  decays  monotonically  for 
large  t,  the  truncation  error  is  an  alternating  series  due  to  the  oscillation 

CO 

of  g(t)exp( jxt) .  Approximation  of  an  alternating  series  J)  (-l)iai,  aj 

i=l  n 

positive,  monotonically  converging  to  0,  by  the  nth  partial  sum  X  (-l)^ai 

i=l 

causes  an  error  bounded  by  aj]4.|.  Choose  T  sufficiently  large  so  that  the 
area  enclosed  by  one  of  the  half-oscillations  is  negligible  compared  to  the 
estimated  integral. 
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We  now  discuss  why  incorporating  Simpson’s  rule  into  the  FFT  gives 
dramatic  improvement  in  our  examples.  This  effect  is  not  obvious,  as  can  be 
seen  in  the  following  case:  if,  in  equation  (1),  f(t)  is  a  highly  oscillatory 
function  such  as  sin(At)  with  A  >>  1,  the  derivatives  increase  in 
magnitude,  and  the  error  bounds  in  Table  1  may  actually  increase  with 
higher-order  integration  methods.  For  example,  if  h  =  1  and  A  >  12  the 
error  bound  in  trapezoidal  rule  is  larger  than  in  rectangular  rule. 

However,  in  our  examples  the  integrals  are  of  the  form  of  equation  (4), 
where  g(t)  is  slowly  varying  in  the  sense  that  its  derivatives  all  have  same 
order  magnitude  (all  bounded  by  M,  say)  on  the  interval  [0,T]. 

Suppose  I(xk)  is  desired  for  equally-spaced  points  Xo  = 
Oi*i»Xaf*iXs  =  First  choose  T  sufficiently  large  so  that  I(.Xk)  is  a  good 
approximation  to  the  original  improper  integral  for  all  x^  in  [0,/?],  and 
each  xk  is  included  in  the  set  {2rTk/T;  k  =  0,1,2,...,N-1}.  For  example, 
suppose  T  =  4n  and  1(0),  1(1),  ...,  1(10)  are  desired.  Our  first  choice  of  N 
is  to  make  N  a  power  of  2  such  that  (N/2)  >  27TkinaLX»  where 
*max  “  2wkoiax/T  =  /?  =  10  (see  the  Nyquist  Criterion  [4]).  In  this  case 
i'max  -  20,  and  N  =  256  becomes  the  first  choice. 

Letting  Ej  denote  the  error  bound  given  by  Table  1  for  rectangular 
rule  applied  to  equation  (4),  we  have 

<  |(TV2)(j(x/N)g(c)  +  (1/N)g'(c))|  (10) 

where  0  <  c  <  T.  For  Simpson’s  rule,  the  error  bound  is 

<  l(T/N)  (T/180)(g^''(c)  +  4jxg~(c)  +  6(jx)"g*(c) 

(11) 

+  4(jx)^g'(c)  +  (jx)‘*g(c))|  . 
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Assuming  that  g  and  the  various  derivatives  of  g  are  all  bounded  by  M, 


<  (T/N)(T/2)M[x  +  1]  (12) 

<  (T/N)*(T/180)Mfl  +  4x  +  6x*  +  4x’  +  x^]  (13) 

Consider  first  the  case  for  which  0  <  x  ^  1.  Then 

E  ThM 
1 

E^  <  Th^Mde/lSO)  <  E^  for  h  =  (T/N)  <  1 
Otherwise,  when  x  >  1 

E^  <  MT(xT/N)  <  MT(l/2)  (16) 


(14) 

(15) 


E^  <  (T/180)(T/N)‘‘m(1  +  4x  +  6x*  +  4x'’  +  x‘] 

<  (MT/180)((6)(5)x‘(T/N)*]  <  MT(l/2)‘/60  <  E^ 


(17) 


since  (xT/N)  <  (1/2). 

For  this  first  choice  of  N  we  therefore  expect  more  accuracy  from 
Simpson’s  rule.  In  actual  practice,  however,  the  choice  of  N  is  determined 
as  is  commonly  done,  by  doubling  the  number  of  partitions  N  (each  step 
reusing  previously  computed  quantities)  until  convergence  is  obtained  in  a 
Rhomberg  scheme.  In  this  process,  each  time  N  is  doubled  the  ratio  (xT/N) 
is  halved,  hence  the  inequality  <  E,  in  (17)  is  strengthened. 

A  single  call  to  FFT  provides  values  for  the  next  iteration  in  Rhomberg’s 
method  (4]  simultaneously  for  all  sample  points  (xj^).  In  this  way,  a 


41 


"parallel-processing"  numerical  integration  algorithm  is  achieved  on  a  serial 
computer  via  the  FFT. 


V.  APPLICATION 

As  an  example  relevant  to  electromagnetics,  suppose  that  the 
Sommerfeld-type  integral 

exp(.ixt)exp(-p|y|  )dt 

2p 

where  p“  =  t*  +  g^,  is  required  at  2M  evenly-spaced  values  xj^  €  [0,2W]. 
Let  T  =  {2nM)L,  where  L  is  an  integer  large  enough  for  adequate 
approximation.  Taking  the  real  part  of  the  integral  allows  replacement  of  the 
lower  limit  by  zero.  By  its  definition, 

Xj^  =  2rTk/T  =  k/LM;  k  =  0,1,2, ...  ,N-1  .  (19) 

Choose  N  for  desired  resolution,  and  recover  every  L’th  evaluation  by  the 
FFT. 

Results  from  this  technique  are  identical  to  those  obtained  by  other 
methods;  moreover,  the  FFT/Simpson's  rule  combination  is  much  faster  since  a 
single  call  to  FFT  returns  the  required  sequence  of  integrals  for  the  set 
(xi(}.  This  example  is  closely  related  to  the  last  example. 

VI.  CONCLUSION 

Numerical  integration  for  aperiodic  functions  with  the  FFT  is  greatly 
improved  by  the  concurrent  use  of  Simpson’s  rule.  Some  algorithms  that 
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already  make  use  of  the  FFT  could  benefit  from  this  slight  modification.  It 
also  appears  to  be  a  promising  technique  for  the  evaluation  of 
Sommerfeld-type  integrals. 
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cols' 


Integration  Rule 


Error 


1.  h  2  f(t  ) 
i=0 


h(b-a)f  •(r7)/2 


-h^(b-a)f  •  •(ti)/12 


Cf(a)  -  f(b)] 


h^(b-a)f^^''^(n)/720 


[f(a)  +  4f(Cj)  -f  2f(t2)  +  .  .  .  +  f(b)] 


-h'*(b-a)f^^'^^{T7)/180 


TABLE  1 

Numerical  integration  methods.  The  rule  names  are: 
trapezoidal,  (3)  corrected  trapezoidal,  and  (4)  Simpson’s. 


(1)  rectangular,  (2) 
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4.3  EM  interactions  among  integrated  microstrip  devices 

This  research  topic  includes  the  scattering  of  microstrip  waves  by  inte¬ 
grated  devices  and  the  large-scale  interactions  among  such  devices.  Inter¬ 
actions  in  the  MMIC  environment  have  been  conceptualized  through  an  EFIE 
description  of  the  microstrip/device  system,  based  upon  the  electric  dyadic 
kernels  exposed  in  Section  4.1.  The  general  electromagnetic  formulation  for 
such  interactions  is  presented  in  [16]. 

Coupling  among  microstrip  devices  can  be  studied  through  a  perturbation 
solution  which  exploits  currents  supported  by  isolated  system  components. 
Significant  scattering  of  microstrip  waves  is  anticipated  only  when  they  in¬ 
teract  with  near-resonant  devices.  The  quasi-closed-form  Chebychev  series 
representations  for  microstrip  line  and  device  currents  detailed  in  Sections 
4. 5-4. 6  permits  the  quantification  of  scattering  coefficients  based  upon  such 
a  perturbation  formulation. 

The  description  of  microstrip-based  circuits  by  integral-operator  equa¬ 
tions  which  exploit  the  electric  Green’s  dyad  kernels  leads  to  a  number  of 
general  properties  and  phenomena  associated  with  such  circuits.  An  excitation 
theory  for  discrete  microstrip  line  and  cavity  modes  is  obtained,  which  yields 
the  excitation  amplitudes  as  overlap  integrals  of  impressed  currents  with  the 
modal  fields  of  those  devices.  A  coupled-mode  perturbation  theory  for  propag¬ 
ation  constant  or  resonant  frequency  shifts  due  to  coupling  among  systems  of 
multiple  microstrip  lines  or  cavities  is  developed,  again  based  upon  the 
integral-operator  description  of  the  sysyem.  Research  paper  [16]  exposes  all 
these  results. 

Large-scale  interactions  among  systems  of  microstrip  devices  can  be 
quantified  by  the  methods  and  approximations  described  in  [16].  The  research 
topics  of  Sections  4. 4-4. 6  are,  in  fact,  specializations  of  such  interactions, 
consequently  [16]  again  provides  the  electromagnetic  foundation  for  the  study 
of  those  topics. 
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Research  Paper  No.  [16] 

DEDUCTION  OF  EM  PHENOMENA  IN  MICROSTRIP  CIRCUITS  FROM  AN  INTEGRAL-OPERATOR 
DESCRIPTION  OF  THE  SYSTEM 

Dennis  Paul  Nyqulst 

Departaent  of  Electrical  Engineering,  Michigan  State  University,  East  Lansing,  Michigan, 
USA  48824 


Ai>3,tr4«;.t 


An  electrlc-fleld  Integral  equation  (EFIE)  description  of  alcrostrlp-based  circuits  permits  the 
relatively  general  deduction  of  EM  phenomena  which  influence  the  behavior  of  those  circuits. 
Prior  efforts  by  many  researchers  have  used  the  EFIE  formulation  to  study  specific  mlcrostrlp 
devices.  The  present  research  exploits  an  EFIE  description  to  study  more  broadly  the  EM  phenom¬ 
ena  and  interactions  In  mlcrostrlp  circuits.  Included  are  the  following  topics;  1]  Identifica¬ 
tion  of  the  complete  propagation-mode  spectrum  for  mlcrostrlp  from  a  singularity  expansion  of 
Its  current.  Including  an  excitation  theory  for  discrete  nodes;  2)  excitation  of  patch-device 
resonances,  which  arise  from  temporal  pole  singularities;  and  3)  formulation  of  a  coupled-mode 
perturbation  theory  leading  to  propagation  constant  and  resonant  frequency  shifts  due  to  coup¬ 
ling  among  mlcrostrlp  lines  and  patch  resonators. 


1.  Introduction 


An  Integral-operator  formulation  Is  presented  for 
electromagnetic  Interactions  In  alcrostrlp-based 
electronic  circuits.  It  Is  a  full-wave  formulation 
based  upon  the  electric  dyadic  Green' s  function  for 
currents  Immersed  In  the  layered  integrated  elec¬ 
tronics  environment.  Full-wave  Integral  equation 
methods  have  been  used  extensively  for  the  study  of 
specific  mlcrostrlp  devices.  A  survey  of  such 
methods  applied  to  the  study  of  discrete  mlcrostrlp 
propagation  modes  Is  Included  In  (II,  while  a  simi¬ 
lar  treatment  for  mlcrostrlp  patch  antennas  Is  ex¬ 
posed  In  (21.  The  present  research  attempts  to 
deduce  more  broadly  those  EM  phenpaena  which  Influ¬ 
ence  the  behavior  of  mlcrostrlp  circuits  from  an 
electrlc-fleld  Integral  equation  description  of  the 
circuit  system. 

2.  EFIE  Description  of  Mlcrostrlp  Circuits 

The  mlcrostrlp  circuit  configuration  consists  of  ' 
conducting  circuit  devices  embedded  In  the  cover 
layer  adjacent  to  the  flla/cover  Interface  of  a 
trl-layered  conductor/flln/cover  environment.  A 
dielectric  film  layer  of  thickness  t  and  refractive 
index  n^.  occupies  -t<y<0  between  a  perfect  conduc¬ 
tor  In  y<-t  and  a  dielectric  cover  of  Index  n  In 

c 

y>0.  Circuit  components,  of  generally  non-vanlsh- 
Ing  thickness,  are  consequently  located  In  yao, 
where  the  y  axis  Is  normal  and  the  x,z  axes  tangen¬ 
tial  to  the  fllm/cover  Interface.  Total  circuit 
surface  S  Is  composed  of  N  components,  where  the 
component  contributed  by  the  n' th  circuit  element 
Is  S^.  Fig.  1  depicts  the  circuit  configuration  In 

cross  section,  and  Indicates  a  general  pair  of  de¬ 
vice  elements. 


OOVEM  LAYEM  n 


CONDUCTOR 


Figure  1.  Configuration  of  mlcrostrlp-based 
electronic  circuit. 


Impressed  current  maintains  Impressed  electric 
field  which  excites  surface  current  ^  on  the 
conducting  circuit  system.  Induced  current  t  main¬ 
tains  a  scattered  field  S*.  The  boundary  condition 


for  tangential  f  at  the  circuit  surface  requires 
t-(6  ♦?*)  «  0.  where  ?  is  a  unit  tangent  vector  at 
any  point  on  that  surface.  Expressing  In  terms 

of  electric  dyadic  Green's  function  C®  for  the  lay¬ 
ered  environment  leads  to  the  electrlc-fleld  inte¬ 
gral  equation  (EFIE)  for  Induced  on  the  circuit 

t-f  5*Crlr' )-^(r' )dS'  «  -  —  g^(r)  (1) 

Js  ^c 

for  all  r«S,  where  k  «n  k„  and  q  «Tt,.''n  with  k^.n,. 

ecu  cOc  00 

the  free-space  wavenumber  a.id  Impedance. 


The  electric  Green’s  dyad  13.4)  can  be  expressed  as 

if(rlr')  -  PV(k^*TO- )5(rtr' 1  ♦  Cs(r-r' )  where  C  Is 
c 

a  Hertzian  potential  Green's  dyad.  Notation  PV 

Indicates  that  should  be  Integrated  in  a  princi¬ 
pal  value  sense  (SI  by  excluding  an  Innate  ‘slice* 
principal  volume  (6)  to  accomodate  the  source-point 

singularity  at  y«y' ,  and  C  >  -yy  Is  an  associated 
depolarizing  dyad.  The  Hertz  potential  Green’s 

dyad  Is  expressed  as  C  •  Tc^  ♦  where  G**  yields 

the  principal  wave  excited  by  ^  in  an  unbounded 
cover  and 


CiFi?') 


A-rA 

xC^x 


,  ac 


-FA 


ac 

r 

dz 


C  A 
Z 


zGtZ 


yields  the  reflected  wave  scattered  by  the. layered 
environment.  G^  yields  potential  components  tang¬ 
ential  to  the  flla/cover  Interface,  while  gJ^  yields 
normal  potential  due  to  normal  current  and  G^  coup¬ 
les  tangential  current  to  normal  potential.  Reci¬ 
procity  of  the  electric  Green’s  dyad  leads  to 
.e 

'«8' 


g:,(?i?')  -  g!  (? 


Pa 


Ir). 


The  several  dyad  components  have  Sommerfeld  inte¬ 
gral  representations 


gP(?i?' )  •  jj  - 


jt. (?-?')  -Pc'y-y'i 


2(2n)^p 


08 

'R^(A) 

■  ff 

R  (A) 

n 

JJ 

n 

,gP(?|?') 

-• 

C(\) 

<?  (?-?') -Pc' , 

5 -  d^x 

2(2it)  p. 
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current  becoaes 


where  t  •  end  d^X  «  d^d<.  Uave- 

nuaber  paraaeters  are  »  (X^-n^^)*^*:  Xe{p^}>0 


for  1  •  f.c.  Reflection  and  coupling  coefficients 
are 

AJX)  A_(X)  2(Nf  -Dp 

Z*(X)  ’  '  2*'(X)2*{X) 

with 

At  •  Pc  -  PfCoth(pj.t)’ 

*n  ■  '^cPc  ■  Pf‘“h(pj.t)' 

Z**  "  Pg  ♦  Pj.coth(pj.t)^ 

2*  -  N^^p^  ♦  Pj.tanh(pj.t) 

where  contrast  factor 

-  "/“c’  Z**’ "  «  0  is 

associated  with  TE  and  TM  surface-wave  poles,  re¬ 
spectively,  of  the  layered  stirround  environaent. 


3.  Propagation-Mode  Spectrua  of  Microstrip 

The  aicrostrip  line  configuration  is  that  speci»i- 
ization  of  Fig.  1  where  the  conducting  eleaents 
have  infinite  extent  along  the  waveguiding  z-axis. 
Infinite  axial  syaaetry  renders  the  Integral  oper¬ 
ator  in  EFIE  (1)  convolutional  in  the  axial  vari¬ 
able,  and  Foiu'ier  transforaation  on  that  variable 
is  consequently  proapted.  Transfom-doaaln  quanti¬ 
ties  are  denoted  by  lower-case  letters  such  that 

9^<^(?,z)>  •  for  all  such  quantities,  where 

p  ■  )W«'ZZ  is  the  2-D  transverse  position  vector  and 
<;  is  the  transfora  variable  corresponding  to  z. 
Subsequent  to  application  of  the  convolution 
theorea,  ^lE  (1)  becoaes 

■•C  ”0 

for  all  ^eC,  where  C  is  the  boundary  contour  about 
the  aicrostrip  conductor  in  the  transverse  plane. 

Green’s  dyad  g*(pl?';<)  ■  PV(k^*S^')g(plp':C)  ♦ 
Cd(p-p' )  where  9  ■  with  the  transverse 

operator  and  g(^|^':<)  >  V^lCCpIp'iz)!.  The  trans- 
fora-doaain  electric  Green' s  dyad  has  the  useful 
reciprocal  property  g^(?l?';G)  •  g^(p'lp:-<). 


i?(?.z)  -  tj[  ^  I  ,  ir(2.0eJ^"d< 

P  ‘"b  (5) 

where  the  sub  of  discrete  modes  Is  contributed  by 
pole  singularities  and  the  continuous-spectrum  con¬ 
tribution  arises  from  deforming  the  Integration 

♦ 

path  about  the  branch  cuts  along  C^.  The  continu¬ 
ous  spectrua  current  k(p, ()  is  consequently  Identi¬ 
fied  as  the  forced  solution  to  EFIE  (2)  for  points 
♦ 

along  CT-  A  discrete-mode  excitation  theory  for 
“  ♦ 

amplitudes  a~  Is  exposed  in  the  next  section. 


■i.  Excitation  and  Coupling  of  Mlcrostrlp  Lines 

The  amplitude  coefficients  a^  and  pole  order  t  In 
representation  (3)  are  established  here.  If  the 

reciprocal  property  of  ^  Is  exploited  subsequent 
to  applying  the  testing  operator  (which  Implements 

the  operation) 


11a  r  dfi^(p)-( - ) 

Wp  Jc  P 


to  EFIE  (2).  then  expansion  of  g*  in  a  Taylor's 
series  about  Is  prompted.  Invoking  defining 

EFIE  (4)  for  k"  leads  to  vanishing  of  the  contribu¬ 


tion  froa  the  first  term  of  that  expansion.  The 
leading  non-vanishing  term  is  consequently  the 
first-order  contribution  proportional  to  (C±C  ). 

^  P 

When  representation  (3)  for  k  Is  used  In  the  re¬ 
sulting  operator  equation.  It  Is  clear  that  the 
leading  term  can  annul  at  most  a  first-order  pole; 
It  is  thus  established  that  (■!. 


If  reciprocity  of  ^  is  invoked  again,  and  only 
leading  terns  are  retained,  propagation-mode  ampli¬ 
tudes  excited  by  e^  are  obtained  as 


where  )^  satisfies  the  normalization  constraint 
P 


Discrete  propagation  nodes  correspond  to  pole  sing¬ 
ularities  of  at  hs*r  which 

■  a^*{?)/((:±<:p)^  (3) 

where  ^p  is  the  eigenaode  current  of  a  wave  propag¬ 
ating  in  the  Iz  direction.  It  is  deaonstrated  in 
Sec.  4  that  f>l,  so  the  poies  are  staple.  If  rep¬ 
resentation  (3)  for  is  exploited  In  ^lE  (2) 

then,  since  ()  is  regular  at  ,  i(^  must  sat- 

P  P 

Isfy  the  hoBogeneous  EFIE 

?(?lp' :<)  •){*(?')-  0  ...  for  all  ^  (4) 

with  non-triviai  solutions  only  for  Th^ 

discrete  propagation  nodes  are  consequently  de¬ 
scribed  by  EFIE  (4). 


-  f  df  f  drf(?)-»l!l|^^  -ilip) 
P  Jc  Jc  P  "  zc  P 


■  ±1 

(7) 


Coupling  coefficient  (6)  Is  an  overlap  Integral  of 
the  iapressed  field  with  the  discrete  eigenmode 

current.  If  e^  maintained  by  7^  Is  expressed  In 

terms  of  g*.  and  the  reciprocity  of  g*  Is  subse¬ 
quently  exploited,  then  an  alternative  coupling 

coefficient  involving  the  overlap  Integral  of 

and  Op  is  obtained.  Using  the  latter  excitation 

amplitudes  In  expression  (S)  leads  to  the  discrete- 
mode  current 


for  (z-z' )«±. 


(8) 


The  discrete  propagation  nodes  corresponding  to 

slnple-pole  singularities  of  i{(p, C)  constitute  one 
contribution  to  the  singularity  expansion  of  that 
current.  Square  root  branch-point  singularities  In 

the  Integral  representation  of  g*  lead  to  slallar 

singularities  in  g*and  ic.  with  associated  branch 
cuts  In  the  complex  C'Plxne.  If  the  real-line  In¬ 
version  integral  is  deformed  Into  the  upper) lower 
half  C-plane  for  z>z'lz<z',  then  the  space-domain 


The  discrete  system  eigenaode  currents  supported  by 
N  coupled  aicrostrip  lines  again  correspond  to 
simple-pole  ^-plane  singularities,  and  satisfy  the 
coupled  system  of  homogeneous  EFIE's 

N 

T  [  g*(plp';<)-ic;;  (p' )df'  »  0  ...  all  p«C  , 

■  hSl  Jc  "P  " 

n  •  for  n«1.2 . M  (9) 

with  non-trivlal  solutions  for  .  For  nearlv- 

P 
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degenerate  coupling,  the  systee-eode  C  Is  near  to 

(01  ^ 
all  the  associated  with  elgennode  currents 

of  the  Individual  Isolated  elcrostrlp.  If 
"P  -e 

reciprocity  of  g  Is  exploited  subsequent  to  appli¬ 
cation  of  the  operator 


<U  - 1 

■P 


for  e-l.Z . N 


to  EFIE's  (9),  then  expansion  of  g  in  a  Taylor's 
series  about  prompted.  The  leading-tern 

contribution  vanishes  for  n^  by  virtue  of  the  de¬ 


fining  EFIE  for 

np 


while  the  first-order  term 

annuls  the  simple  poles  of  i^(p,  C).  Retaining  only 
the  leading  non-vanishing  terms,  subsequent  to  the 


perturbation  approximation  k' 


_  .  leads  to 

np  n  np 

the  homogeneous  system  of  algebraic  equations 


pole,  which  establishes  f«l.  If  representation 
(12)  Is  used  in  the  resulting  expression,  and  reci¬ 
procity  is  again  invoked,  then  the  excitation  amp¬ 
litude  is  obtained  as  the  overlap  Integral 

Jk  r  1 

a  (u)  - - ^  r  iJ  (r)-r(r,u)dS  (14) 

P  "cS  P 

where  normalization  constant  C  Is 


C 

P 


aflrlr*  :u) 
do 


■t  (r')dS'.  (15) 
"p  P 


An  EFIE-based  coupled-mode  perturbation  theory  for 
a  nearly-degenerate  coupled  system  of  N  resonant 
devices  is  based  upon  the  coupled  system  (13)  of 
IE's.  A  procedure  analogous  to  that  which  lead  to 
(10)  above  yields  the  algebraic  equations 

Iu-«‘°’]a  *Tc  a  -0...  for  n«l,2 . N  (16) 

p  fli  L  An  n 
^  n»A 


♦  r  C^a;  a  0  ...  m-1,2 . N  (10) 

with  non-trivtal  solution  only  when  det«l  >  0. 
Coupling  coefficients  are  given  by  overlap  inte¬ 
grals  of  the  m' th  strip  current  with  the  n'th  strip 
field  as 


C 


i 


mn 


g»(0)(^).dt(0)(d,^^ 

Ap  np 


(11) 


while  normalization  coefficients  are  given  by 

AA 

expression  (7)  applied  to  Isolated  eigenmode  cur¬ 
rents  fiO) 

on  the  m' th  strip. 


For  nearly  degenerate  coupling  of  two  microstrip, 
the  perturbation  approximation  to  system-mode  prop¬ 
agation  eigenvalues  leads  to  <  >  <  > 

is  the  average  of  isolated  phase  con¬ 
stants  and  4  ■  (d^+K^)*^  whore  d  «  (<j-<2)/2  and 


5.  Excitation  and  Coupling  of  Patch  Resonators 

Currents  excited  on  a  resonant  patch  device  by  im¬ 
pressed  radiation  satisfy  EFIE  (1),  where  S  Is  de¬ 
scriptive  of  a  single  patch  or  a  coupled  system  of 
such  resonant  devices.  Resonant  currents  are  as¬ 
sociated  with  pole  singularities  at 

temporal  frequency  domain,  such  that 

i!(?,u)  a  a  (u)i!  (?)/(«-«  )^  .  . .  o  near  u  (12) 

P  P  P  P 

where  it  (?)  is  the  normalized  resonant  elgencur- 
P 

rent.  Subsequent  to  using  representation  (12)  in 

EFIE  (1),  and  recognizing  that  is  regular  at  o  , 
^  P 

it  is  evident  that  must  satisfy  the  homogeneous 

EFIE 

t  f  C®(?l?':o)  it  (?')dS'  -  0  ...  all  reS  (13) 
■'S  ** 

with  non-trlvlal  solution  nly  for  o^^.  Equation 
(13)  defines  the  resonant  patch  eigenmode. 

If  reciprocity  of  C*  is  exploited  subsequent  to 
application  of  the  operator  dS  i{p(r')'( - )  to 

EFIE  (1).  then  expansion  of  &  in  a  Taylor's  series 

about  u~u  is  prompted.  The  first  term  of  the  re- 
P 

suiting  expression  vanishes  by  virtue  of  definition 
(13)  for  iJp,  so  the  leading  hon-vanlshlng  term  Is 

linear  In  (w-w  );  it  can  annul  at  most  a  simple 

P  • 


which  yields  the  system-mode  resonance  frequency  as 

the  solution  to  det(u]  -  0,  where  are  the  iso- 

P 

lated  resonant  frequencies.  Coupling  coefficients 
are  given  by  the  overlap  Integral  of  Isolated  cur¬ 
rents  and  fields 

C  -^f  ii^°’(r)-g'^°’(r)dS  (17) 

■n  Jg  ap  np 

a 

where  currents  are  normalized  by  equating  coeffic¬ 
ients  similar  to  (IS)  to  unity. 

6.  Concluding  Remarks 

Several  general  properties  of  Bi  phenomena  in  mi- 
crostrlp-based  circuits  have  been  deduced  from  an 
Integral-operator  description  of  the  circuit  sys¬ 
tem.  The  complete  propagation-mode  spectrvin  of  a 
microstrip  line  was  Identified  from  a  singularity 
expansion  of  its  current.  It  is  believed  tliat  this 
is  the  first  such  treaUwnt  of  the  continuous  spec¬ 
trum.  An  excitation  theory  for  discrete  nodes 
along  Isolated  microstrip  and  a  coupled-mode  per¬ 
turbation  theory  for  nearly-degenerate  coupling 
among  any  number  of  microstrip  were  exposed.  Patch 
device  resonances  were  identified  with  temporal 
pole  singularities,  while  excitation  and  perturba¬ 
tion  coupling  theories  for  such  resonant  elements 
were  advanced.  Typical  niiaerlcal  results  will  be 
included  in  the  oral  exposition  of  the  paper. 
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4.4  Propagation-mode  spectrum  of  microstrip  guiding  structures 

Research  on  the  complete  propagation-mode  spectrum  of  microstrip  guiding 
structures  exploits  a  transform-domain  integral-operator  description  of  the 
structure  and  is  based  upon  a  singularity  expansion  of  the  microstrip  current 
[16].  It  includes  an  excitation  theory  for  discrete  and  continuous  spectral 
components.  Discrete  propagation  modes  correspond  to  pole  singularities  in 
the  complex  transform  plane,  while  the  continuous  spectrum  arises  from  inte¬ 
gration  about  branch  cuts  associated  with  branch-point  singularities  of  the 
Green’s  dyad  kernels  [2,5,7,16]. 

The  discrete-mode  quantification  has  been  generalized  [13,15]  to  include 
the  possibility  of  a  ferrimagnetic  thin-film  layer.  This  study  is  based  upon 
a  transform-domain  EFIE  description  for  current  induced  on  the  microstrip;  the 
dyadic  Green’s  function  kernel  described  above  is  exploited.  The  discrete 
(principal,  and  several  higher-order)  modes  are  now  quite  well  understood.  A 
new  Chebychev-Galerkin’ s  MoM  numerical  solution  has  been  implemented  using 
Chebychev  polynomials,  modulated  by  appropriate  square-root  edge-condition 
factors,  for  expansion  and  testing  of  the  transverse  dependence  of  unknown 
eigenmode  currents  (both  longitudinal  and  transverse). 

Results  indicate  that  the  transverse  current  component  is  negligible  for 
the  principal  discrete  eigenmode,  but  becomes  important  in  the  higher-order 
modes.  This  method  provides  a  rapidly  convergent  quasi-closed-form  solution 
for  the  eigenmode  currents,  since  analytical  expressions  for  the  current  com¬ 
ponents  are  available  subsequent  to  numerical  quantification  of  expansion 
coefficients  in  the  Chebychev  series.  Convergence  is  so  rapid  that  two  terms 
suffice  for  the  principal  mode,  and  at  most  four  need  be  retained  for  even  the 
second  higher-order  mode.  It  was  found  that  coupling  to  surface-wave  modes  of 
the  layered  surround,  and  consequent  leakage  phenomena,  is  unlikely  to  occur 
for  practical  configurations  with  dielectric  thin-film  layers.  Presence  of  a 
ferrimagnetic  film  layer  may  alter  this  conclusion,  however,  and  numerical 
studies  remain  to  be  completed.  The  quasi-closed-form  solutions  for  discrete 
eigenmode  currents  are  particularly  suitable  for  application  in  perturbation- 
type  studies  of  their  coupling  with  adjacent  devices  and  other  nearby  micro¬ 
strip  lines.  Extensive  numerical  results  are  presented  in  connection  with  the 
formulation  for  coupled  microstrip  lines  exposed  in  Section  4.5. 

Radiation  fields  of  the  continuous  spectrum  have  been  conceptually  iden- 
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tified  as  forced  solutions  to  the  EFIE’ s  at  points  along  transform-plane 
branch  cuts.  Their  numerical  quantification  has  been  ititiated  through  a  MoM 
formulation  [7]  for  the  forced  microstrip  currents.  Additional  research  on 
potentially  efficient  approximate  iterative  solutions  remains  to  be  completed. 
The  radiation  field  is  required  for  any  serious  investigation  of  scattering  by 
discontinuities  along  open  microstrip  in  MMIC’ s.  Since  no  suitable  represen¬ 
tation  of  the  radiation  field  in  microstrip-based  circuits  is  presently  avail¬ 
able,  this  topic  constitutes  one  of  the  major  unsolved  problems  in  the  theory 
of  such  circuits.  It  is  conjectured  that  a  suitable  formulation  can  be  ob¬ 
tained  based  upon  the  integral-operator  description  of  microstrip,  since  that 
full-wave  theory  at  least  provides  for  a  conceptual  definition  of  the  radi¬ 
ation  field  and  its  spectrum  [7,161. 

Research  papers  [2,5,7]  follow  on  the  subsequent  pages,  while  papers 
[13,15]  were  included  in  Section  4.1  and  paper  [16]  was  appended  to  Section 
4.  3. 
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Research  Paper  No.  [2] 


RADIATIVE  AND  SURFACE  WAVE  LOSSES  IN  MICROSTRIP  TRANS¬ 
MISSION  LINES:  J.  S.  Bagby,  Department  of  Electrical 
Engineering,  University  of  Texas  at  Arlington, 
Arlington,  Texas  76019,  and  0.  P.  Nyquist,  Department 
of  Electrical  Engineering  and  Systems  Science,  Michigan 
State  University,  East  Lansing,  Michigan  48824. 

An  exact  dyadic  Integral  equation  Is  utilized  In  the 
analysis  of  propagation  In  uniform  Integrated  mlcrostrlp  trans¬ 
mission  lines.  The  object  of  the  analysis  Is  to  predict  and 
quantify  the  radiative  and  surface  wave  losses  In  such  systems. 

The  axial  ly-trans formed  surface  current  i<(p;C)  of  a 
natural  mode  on  an  Integrated  mlcrostrlp  line  satisfies  the 
honogeneous  Integral  equation: 

t«(k^  +  vv*)  ^g(?|p';c)*i<(p';C)  dt'  =  0,  pci 

where  k  Is  the  wavenumber  in  the  cover  region,  g(p!p';c)  Is 
the  Hertzian  potential  Green's  dyad  of  the  background  struc¬ 
ture,  c  Is  the  complex  propagation  constant  of  the  mode,  t  Is  a 
unit  tangent  to  the  transmission  line,  and  l  Is  the  cross- 
sectional  contour  of  the  line. 

The  Integral  representation  of  the  Green's  dyad  exhibits 
singularity  when  the  real  part  of  c  Is  less  than  the  propag¬ 
ation  constant  of  a  surface  wave  mode  In  the  Integrated  circuit 
background  structure.  It  Is  shown  that  for  the  dominant  mode 
of  the  line  the  singularity  is  not  Implicated,  and  all  losses 
are  due  to  radiation  Into  the  cover  medium.  For  higher  order 
modes,  however,  the  singularity  Is  Implicated,  requiring  Inclu¬ 
sion  of  a  residue  term  representing  excitation  of  surface  waves 
In  the  film  layer.  For  higher  order  modes,  higher  losses  are 
observed  due  to  excitation  of  surface  waves  In  the  film  layer. 

This  equation  Is  solved  numerically  by  the  method  of 
moments  for  narrow  and  wide  mlcrostrlp  transmission  lines.  The 
above-described  effects  are  demonstrated  In  both  cases. 


cover  medium,  Gc 
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Research  Paper  No.  [5] 


COMPLETE  PROPAGATION-MODE  SPECTRUM  OP 
MICROSTRIP  GUIDING  STRUCTURES 

M.J.  Cloud  and  D.P.  Nyquist 
Department  of  Electrical  Engineering 
and  Systems  Science 

B.C.  Drachman,  Department  of  Mathematics 
Michigan  State  University 
East  Lansing,  Michigan  48824 


Results  of  investigating  the  propagation-mode 
spectrum  for  longitudinally-invariant  microstrip 
systems  are  presented.  The  influence  of  a  layered 
conductor /film/ cover  background  environment  is 
rigorously  accounted  for. 

Surface  current  eigenmodes  K(r)  are  solutions  of 
a  homogeneous  electric  field  integral  equation 
(EFIE).  The  EFIE,  which  is  enforced  at  all  points  r 
on  the  conducting  microstrip  surface  S,  may  be  stated 
as 


^  .  (k^  +  V  7  •)  /  '5(r,r'  )  •  K(r'  )  dS*  =0 

J  S 

where  k  is  the  wavenumber  in  the^c§ver  medium,  ^  is  a 
unit  tangent  vector  to  S,  and  G(r,r')  is  a  Hertz 
potential  Green's  dyad.  This  dyad  depends  explicitly 
upon  the  thickness  and  permittivity  of  the  film 
layer.  Often  neglected  transversely-flowing  surface 
currents  are  retained  in  the  model  through  a  coupled 
integral  equation  approach,  and  are  quantified  along 
with  the  longitudinal  components. 

Emphasis  is  placed  on  leaky  microstrip  modes 
associated  with  transversely  propagating  surface 
waves  of  the  background  environment.  Numerical 
soldtion  is  performed  by  Galerkin's  method  using 
entire-domain  basis  functions. 
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Research  Paper  No.  [7] 


Ph.D.  DISSERTATION  ABSTRACT 


by 


Michael  John  Cloud 

Department  of  Electrical  Engineering 
Michigan  State  University 
East  Lansing,  Michigan  48824 

December  1987 


Electromagnetic  waves  in  millimeter-wave  integrated  circuits  are  studied. 
Dielectric  film  and  cover  regions  overlay  a  conducting  half-space  in  the  con¬ 
figuration  modeled,  forming  a  non-uniform  background  for  integrated  devices. 
Emphasis  is  placed  on  the  microstrip  transmission  line  in  this  environment. 

This  research  exploits  an  integral-operator  description  of  the  system. 
Constructed  through  Hertz  potentials,  an  electric  field  integral  equation  (IE) 
quantifies  microstrip  surface  currents.  Complex  plane  analysis  in  the  axial 
Fourier  transform  domain  leads  to  the  rigorous  identification  of  discrete  and 
continuous  eigenvalue  spectra.  Discrete  modes  are  associated  with  simple  pole 
singularities,  while  the  continuous  spectrum  arises  from  branch-cut  integrals. 
These  modal  spectra  are  subsequently  linked  to  natural  and  forced  solutions  of 
the  IE. 

Discrete  wave  modes  are  associated  with  the  homogeneous  IE,  which  is 
solved  by  the  moment  method  for  eiectrically-thin  microstrip.  The  current 
density  function  is  expanded  in  both  subsectional  and  entire-domain  basis 
functions.  Results  clearly  validate:  (1)  the  dominant  axial  current  approxi¬ 
mation  for  narrow  strips,  and  (2)  the  edge  singularity  for  axial  currents. 
Dispersion  characteristics  are  given  for  several  modes;  those  for  the  fundam¬ 
ental  mode  compare  favorably  with  results  already  in  the  literature. 

Solution  of  the  forced  IE,  at  points  along  the  complex-plane  branch  cut, 
yields  spectral  components  of  the  strip  radiation  modes.  The  forcing  function 
is  taken  to  be  the  electric  field  impressed  by  a  vertical  monopole  current 
that  resides  in  the  film  layer.  Preliminary  moment-method  results  are  given. 

The  transform-domain  IE  is  also  a  basis  for  the  study  of  coupling  phenom¬ 
ena  in  axially-uniform  multi-strip  systems.  System  propagation  modes  are 
characterized  by  coupled  currents  sharing  simple-pole  singularities.  Pertur¬ 
bation  approximations  apply  for  loose  coupling,  leading  to  an  overlap-integral 
description  of  the  coupling  phenomenon. 

Finally,  a  novel  approach  to  the  numerical  evaluation  of  Sommerfeld  inte¬ 
grals,  using  the  Fast  Fourier  Transform,  is  advanced. 
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4.5  Coupling  between  adjacent  microstrip  lines 

Discrete  system  modes  of  a  coupled  multi-microstrip  configuration  can  be 
studied  based  upon  a  transform-domain  integral-operator  description  of  the 
system.  Those  system  modes  arise  fron  pole  singularities  in  the  complex 
transform  plane.  A  Chebychev-Galerkin’ s  solution  to  the  homogeneous  coupled 
EFIE’s  for  discrete  eigenmode  currents  of  an  adjacent  microstrip  system  has 
been  developed  [10].  A  manuscript  [18]  describing  this  accurate  and  rapidly- 
convergent  full-wave  formulation  is  presently  in  the  review  cycle.  The  iso¬ 
lated-strip  eigenmode  propagation  constant  is  found  to  split  when  two  adjacent 
microstrip  become  coupled,  leading  to  two  new  system  modes  having  symmetric 
and  antisymmetric  parity  among  the  associated  currents  on  the  two  lines.  This 
formulation  again  leads  to  concise  quasi-closed-form  solutions  for  the  eigen¬ 
mode  currents  supported  by  coupled  microstrip.  Included  in  [18]  are  extensive 
numerical  results  for  both  isolated  and  coupled  microstrip  configurations. 

An  EFIE-based  coupled-mode  perturbation  theory  for  coupling  among  any 
number  of  adjacent  microstrip  lines  has  been  formulated  [11,14,16].  It  leads 
to  a  simple  formulation  for  propagation  eigenvale  shifts  and  relative  ampli¬ 
tudes  of  the  various  line  currents  for  discrete  system  modes.  Since  the  gen¬ 
eral  Green’s  dyad  kernels  are  exploited,  both  dielectric  and  ferrimagnetic 
thin-film  layers  can  be  accomodated.  The  simple  quasi-closed-form  Chebyshev 
series  expressions  for  the  isolated  eigenmode  currents  are  well  suited  for 
application  in  this  coupled-mode  theory.  Numerical  implementation  leads  to 
results  which  agree  almost  exactly  with  those  of  the  numerical  MoM  solution 
except  when  the  lines  are  very  closely  spaced.  A  research  paper  manuscript 
describing  this  integral-operator-based  perturbation  formulation  is  presently 
in  preparation. 

Papers  [10,11,14,18]  are  appended  here,  while  paper  [16]  was  included  in 
Section  4. 3. 
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Research  Paper  No.  [10] 


ENTIRE-BASIS  MOM  ANALYSIS  OF  COUPLED  MICROSTRIP 

TRANSMISSION  LINES 

C.-H.  Lee,  J.S.  Bagby 

Department  of  Electrical  Engineering 

University  of  Texas  at  Arlington 

Arlington,  Texas  76019 

Y.  Yuan,  D.P.  Nyquist 

Department  of  Electrical  Engineering  and  System 
Science 

Michigan  State  University 
East  Lansing,  Michigan  48824 

A  system  of  exact  dyadic  integral  equations  is  utilized  in  the 
analysis  of  propagation  in  uniform  coupled  integrated  microstrip 
transmission  lines.  Axially-transformed  natural  mode  surface 
currents  on  N  coupled  microstrip  transmission  lines  satisfy  the 
homogeneous  coupled  dyadic  integral  equations: 

N  /*..«-  —  - 

I  g  (pip*;0*ki(p*)de*=o,  j=i . N 


where  kj  is  the  transformed  surface  current  on  the  i*'’  strip,  9  i  s 
the  electric  Green's  dyad  of  the  background  structure,  ^  is  the 

unknown  propagation  constant  of  the  coupled  mode,  tj  is  a  unit 
tangent  to  the  strip,  and  t-,  is  L'O  cross-sectional  contour  of  the 
i*'’  strip. 

The  above  integral  equation  is  solved  by  the  method  of 
moments.  Longitudinal  and  transverse  currents  are  expanded  in 
entire  basis  functions  consisting  of  Tchebyshev  polynomials  of  the 
first  and  second  kind  with  multiplicative  factors  incorporating  the 
proper  edge  behavior.  This  formulation  is  shown  to  converge  to 
accurate  results  in  as  few  as  three  terms. 

Numerical  results  in  the  form  of  dispersion  curves  and 
surface  current  distributions  are  presented  for  the  dominant  and 
first  three  higher  order  modes,  both  even  and  odd,  for  the  case  of 
two  identical  coupled  lines. 
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Research  Paper  No.  [11] 


EFIE-BASED  PERTURBATION  APPROXIMATION 
FOR  COUPLED  MICROSTRIP  LINES 
Yi  Yuan  and  Dennis  P.  Nyquist 
Department  of  Electrical  Engineering 
Michigan  State  University 
East  Lansing,  MI  48824 

To  solve  the  eigenvalue  problem  of  coupled  microstrip 
lines,  a  full-wave  analysis  based  on  an  electric  dyadic 
Green's  function  is  developed.  The  electric  field  integral 
equations  (EFIE's)  are  solved,  by  a  Galerkin's  MoM  technique 
with  Chebychev  polynomial  basis  functions,  for  both  isolated 
and  coupled  microstrip.  The  direct  numerical  solution  be¬ 
comes  very  time  consuming  for  a  coupled  microstrip  system, 
consequently  an  approximate  coupled-mode  perturbation  formu¬ 
lation  is  pursued. 


In  this  paper,  we  present  an  EFIE-based  perturbation  ap¬ 
proximation  to  solve  for  the  system  eigenmodes  of  N  coupled 
microstrip  lines.  In  solving  the  coupled  EFIE's,  the  eigen¬ 
mode  current  of  the  isolated  line  (which  is  obtained  by  a 
Galerkin's  MoM  solution  to  the  isolated  EFIE  in  a  convenient 
Chebychev  polynomial  series)  is  used  as  a  first-order  pertur¬ 
bation  approximation  for  nearly-degenerate  eigenmode  currents 
of  the  loosely-coupled  system.  The  EFIE's  yield  a  matrix 
equation 


C  a 

mm  m*’^  ^mp 


C  a 
mn  n 


»  0 


n;£m 


for  msl,2, . . . ,N 


where  C  is  the  unknown  propagation  eigenvalue,  is  that 

Qlp 

eigenvalue  for  the  p'th  mode  on  the  n'th  isolated  microstrip, 
a  is  the  current  amplitude,  and  the  C  are  coupling  coef- 

Q  Bui 

ficients  involving  field/current  overlap  integrals.  C  is 

that  value  which  leads  to  a  non-trivial  solution  for  the  a  . 

n 

Since  the  C  are  C~independent ,  the  numerical  procedure  is 
mn 

relatively  efficient.  For  the  special  case  of  two-line  coup¬ 
ling,  it  is  found  that  the  propagation  eigenvalues  split  and 
shift  symmetrically  away  from  their  isolated  limit  as  the  two 
microstrip  become  closely  spaced. 


A  numerical  implementation  of  the  perturbation  method  is 
developed.  The  numerical  results  of  the  perturbation  approx¬ 
imation  are  compared  with  those  of  the  HoM  numerical  solution. 
Computation  times  are  also  compared,  and  the  validity  range 
of  the  perturbation  approximation  is  investigated. 
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Research  Paper  No.  [14] 


COUPLED  MICROSTRIP  TRANSMISSION  LINES:  FULL-WAVE 
PERTURBATION  THEORY  AND  EXPERIMENTAL  VALIDATION 
* 

Yi  Yuan  ,  John  Vezmar,  Gregory  King  and  Dennis  Nyquist 
Department  of  Electrical  Engineering 
Michigan  State  University 
East  Lansing,  Michigan  48824 


The  coupling  between  adjacent,  parallel  microstrip  transmis¬ 
sion  lines  in  the  micro/mm-wave  PC/IC  environment  is  studied  both 
analytically  and  experimentally.  A  numerically  efficient  coupled- 
mode  theory  is  obtained  through  a  full-wave,  EFIE-based  perturba¬ 
tion  approximation.  The  perturbation  theory  leads  to  propagation 
eigenvalue  shifts  associated  with  sirnmetric  and  antisymmetric 
coupling  of  both  the  principal  and  higher-order  propagation  nodes. 
Experimental  measurements  are  made  on  PC -board  implementations  of 
coupled  microstrip  pairs  at  nicro/mm  wavelengths  to  confirm  the 
predicted  coupling  lengths. 


A  full-wave  EFIE  formulation  for  the  currents  on  N  coupled 
nicrostrip  lines,  located  in  the  cover  layer  at  the  film/cover 
interface  of  a  tri-layered  conductor/film/ cover  environment,  is 
based  upon  the  Sommerfeld-integral  representation  of  an  appropri¬ 
ate  electric  Green's  dyadic.  The  perturbation  solution  to  those 
coupled  EFIE's  exploits  a  rapidly-convergent ,  modulated  Chebychev 
series  numerical  solution  for  the  isolated  strip  currents  as  the 
zero ' th-order  approximate  induced  current  in  the  coupled  system. 
Propagation  eigenvalues  C  Cor  coupled  system  modes  satisfy 


C  a  + 

mm  m  mp 


C  a  =  0 
mn  n 


. . .  for  m  =  1,2, ... ,N 


n^ 


.(0) 


where  C""  is  the  phase  constant  for  the  p'th  mode  on  the  m'th 
mp 

isolated  microstrip,  a  is  the  current  amplitude  on  the  n'th 

strip,  C  is  a  normalization  constant,  and  the  C  are  coupling 
mm  mn 

coefficients  involving  field/current  overlap  integrals.  For  the 
case  of  two  coupled  lines,  the  propagation  eigenvalues  split  and 
shift  symmetrically  away  from  their  isolated  limit  as  the  micro¬ 
strip  become  closely  spaced.  A  simple  coupled-mode  theory  leads 
to  the  coupling  length  in  terms  of  the  eigenvalue  shifts. 


The  coupling  length  is  measured  experimentally,  at  micro/mm 
wavelengths  on  PC-board  implementations  of  two  coupled  microstrip, 
to  validate  the  perturbation-theory  predictions  for  similar  sys¬ 
tem  parameters.  A  monopole  microcoax  probe  is  used  to  measure  the 
axial  electric  field  amplitude  distribution  along  the  coupled  sys¬ 
tem.  The  probe  is  inserted  through  the  conducting  ground  screen 
into  the  film  and/or  cover  layers  to  measure  the  normal  electric 
field  there.  A  pattern  of  small  holes  is  drilled  through  the  PC 
board  to  implement  probe  insertion.  A  harmonic  heterodyne  detec¬ 
tor  assures  measurement  linearity.  A  new  theoretical  study  of 
higher-order  mode  coupling  is  included. 
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ENTIRE-DOMAIN  BASIS  MOM  ANALYSIS  OF 
COUPLED  MICROSTRIP  TRANSMISSION  LINES 


C.-H.  Lee,  J.  S.  Bagby 
Department  of  Electrical  Engineering 
University  of  Texas  at  Arlington 
Arlington,  TX  76019 


Y.  Yuan,  D.  P.  Nyquist 
Department  of  Electrical  Engineering 
Michigan  State  University 
East  Lansing,  MI  48824 


Abstract 


A  full-wave  spectral-domain  integral  equation  formulation  is  used  to 
analyze  coupled  microstrip  transmission  lines.  A  method  of  moments 
solution  is  implemented  utilizing  entire-domain  basis  functions  which 
incorporate  appropriate  edge  conditions  for  transverse  and  longitudinal 
current  components,  allowing  for  closed-form  evaluation  of  spatial  integrals. 
In  contrast  with  earlier  subdomain  basis  solutions,  improved  accuracy  is 
obtained  using  far  fewer  terms.  Numerical  results  in  the  form  of 
propagation  contants  and  current  distributions  are  presented  for  the 
dominant  and  first  two  higher-order  modes,  and  are  compared  to  results  of 
other  techniques. 
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1.  Introduction 


The  problem  of  analysis  of  coupled  microstrip  transmission  lines  is  one  of 
both  considerable  practical  interest,  and  of  long  history.  Traditionally, 
various  quasi-static  methods  [1-5]  have  been  used  to  compute  the  propagation 
characteristics  of  the  coupled  system.  However,  such  methods  are  inherently 
inaccurate  at  higher  frequencies,  and  are  also  found  to  be  inadequate  at  low 
frequencies  for  many  useful  combinations  of  substrate  thickness  and 
dielectric  constant  [6].  In  such  cases  a  more  accurate  full- wave  analysis  must 
be  utilized. 

Recently,  an  exact  spectral-domain  coupled  integral  equation  formulation 
has  been  introduced  [8].  A  subdomain  basis  method  of  moments  solution  of 
the  coupled  equations  has  been  implemented  and  yields  accurate  results  in  the 
form  of  current  distributions  and  propagation  constants  for  the  dominant  and 
several  higher-order  modes  [7].  However,  this  method  requires  a  large 
number  of  basis  functions  to  achieve  desired  accuracy,  yielding  long 
computation  time. 

In  this  paper  an  entire-domain  basis  method  of  moments  solution  is 
utilized.  The  basis  functions  are  chosen  carefully  to  incorporate  the  proper 
behavior  at  the  edges  of  the  strips.  Besides  allowing  closed-form  evaluation 
of  spatial  integrals,  it  is  found  that  good  accuracy  is  obtained  (in  the  case  of 
an  isolated  microstrip)  for  as  few  as  three  basis  functions  for  each  current 
component  of  low-order  modes,  yielding  a  compact  approximate  closed- 
form  solution  for  currents  and  a  substantial  savings  in  computation  time. 
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In  section  2  the  mathematical  formulation  and  numerical  solution 
technique  are  introduced.  Numerical  results  in  the  form  of  propagation 
constants  and  strip  current  distributions  are  presented  in  Section  3. 
Conclusions  and  areas  for  further  investigation  are  given  in  Section  4. 
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2.  Formulation 


Consider  the  coupled  microstrip  geometry  depicted  in  Figure  1.  The 
integral  equation  satisfied  by  the  unknown  surface  current  with  an  assumed 
propagation  dependence  of  exp[j(o)t-Cz)]  on  an  isolated  microstrip 
transmission  line  is  given  in  [8].  This  equation  can  easily  be  extended  to  an  N 
coupled  microstrip  transmission  line  structure  [7],  and  takes  the  following 
form: 


ti<(  V  V-  +k 


(plp';C)-ki(p')de’  =  0,  pee^.  j=l . N. 


(1) 


where  iq  is  the  axially-transformed  surface  current  on  the  i^^  strip,  g  is  the 
transformed  electric  Hertzian  potential  Green’s  dyad  of  the  microstrip 
background  structure,  ^  is  the  unknown  propagation  constant  of  the  coupled 
mode,  *tj  is  a  unit  tangent  to  the  jdi  strip,  t.  is  the  cross-sectional  contour  of  the 

i*  strip,  and  V  =  Vj+  z=x —  +  y —  +  jCz  is  the  transformed  del  operator. 

dx  dy 

The  Hertzian  potential  Green's  dyad  decomposes  into  a  principal  and  a 
reflected  part,  g  =  fg‘’+g^ ,  where  9  is  the  2-dimensional  unbounded-space 
Green's  function  in  integral  representation,  and 

g  (pip')=xg,x+y(— 

dx 

The  scalar  components  of  the  reflected  Green's  dyad  are  given  in  terms  of 
inverse  transform  integrals: 


9c^+gny+jCgcZ)+zg,z  (2) 
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(3) 


222  2^2  2^2 

Here  X  =5  ,Pc  =X  -k,.  ,Pf=X  -kf,  and  the  reflection  and  coupling 
coefficients  are  given  by 


Pc-pfCOth(pft)  ^  Kpc-pftanh(pft) 
‘pc+pfcoth(pf0  ""Kpe+pftanh(pft) 

_ 2(K>l)Pc _ 

[Pc+pfCOth(pft)][Kpc+pftanh(pfi)]  ’ 


(4) 


2.1  Two  identical  microstrip  lines 

Now  consider  a  coupled  system  of  two  identical  microstrip  lines.  Assume 
that  the  lines  are  infinitely  thin  and  are  extended  axially  to  infinity,  as  shown 
in  Figure  2,  where  b  is  half  of  the  spacing  between  the  centers  of  the  two 
strips  and  w  is  half  of  the  strip  width. 

In  this  case,  the  integral  equation  becomes 
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V  y™Z^jg(k/+VV  )g(plx',y'=0;O-k;j(x',y'=0;C)<ix'  =  0,  j=l,2 


i.U‘ 


where  interchange  of  differentiation  and  integration  is  valid  since  y  0.  As 
y->0,  the  contour  involved  in  the  integration  is  only  the  path  along  the  x-axis 
and  is  denoted  as  ^-x..  For  the  specialized  structure,  the  tangential  unit  vector 

t  and  the  surface  current  k  can  be  expressed  as 


t=xtx+ztz,  k=xkx+zk2 


To  evaluate  equation  (5),  we  rewrite  the  integral  representations  of  the  scalar 


components  g**  and  g[,n.<-as 


g'n.o^P  I  =  J  i'n.ctp  I 


with 


J§(*-x').-Pclyl 

.pc  C 

- 


j^(x.x-)  -p^y 


where  the  fact  y'=0  has  been  applied.  Substitute  (6)  into  (5)  to  obtain 


(kc^+V  V  • )  i (p  1 X'; ; k  (x';  Q  dx*  =  0 


(7) 


^  »  9  ^  d 

Notice  that  the  operator  V  is  now  V  =j^xH — y+j^z  since - when  we 

dy  dx 

interchange and  d^. 


After  some  algebraic  manipulation,  we  obtain 


t-  i-  k=(i'’+i[)(t.k,+t,k,) 


.P  .r 


t-  vv  •i-k=(j^t,+jCt,)[(r+i;)(j^,+jCk,)-pX(i^k,+ji;k,)] 


(8) 


(9) 


Substituting  (8)  and  (9)  into  (7)  and  taking  the  limit  y->  0  in  the  results,  we 
obtain  two  coupled  equations  by  letting  t=x  and  t=z.  They  are 
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[k,^(l+Rt)  k^-  ^(1  +  Rt-  PeC )  (^  k^+ c  kj]  dx'  =  0 


[ke\l+Rt)  k^.  C(1  +  Rf  PcC )  (^  k^ + c  k,)]  d^  dx'  =  0 


The  final  form  of  the  integral  equations  for  the  two  indentical  microstrip 
line  structure  can  be  obtained  by  resolving  as  [  -  b  -  w  ,  -  b  +  w]  and 
.  b  +  w] , 


[k,^Rk,i-^(R-C’)(^k,i+Ck,i)]ci^dx' 


(10) 
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-b  +w 


z: 


[k/Rk,i-C(R-C*)(^k,i+Ck^i)]dUx' 


-b  -  w 


b+w 


[k,^R  k,2-  C(R  -  C’)  ( ^  kx2+  C  k^2)]  d^  dx'=  0 


(11) 


Here  kxl,  k2i  and  kx2»  kz2  are  transverse  and  axial  currents  for  strip  1  and 
strip  2,  respectively,  and 


R= 


1  +  Rt= 


2Pc 

p<.+  PfCoth(pft) 


C'  =  PcC  = 


2(K-l)p, 


[Pc+PfCoth(pft)]  [Kpc+Pftanh(pft)] 


Equations  (10)  and  (11)  are  the  coupled  integral  equations  for  the 
microstrip  structure  depicted  in  figure  2.  They  have  a  non-trivial  solution 
for  VijJs  and  k^’s  only  when  C  =  C  m  is  an  eigenvalue;  here  m=0,l,2,...,  is  the 
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mode  number.  In  practical  analysis,  we  can  take  advantage  of  the  symmetry 
in  X  inherent  in  the  coupled  stmcture  to  further  simplify  the  I.E.’s.  For 
symmetric  eigenmodes  the  currents  exhibit  even  symmetry  and  we  let 

kxi(x’)  =  kx2(-x')  andk2i(x')=  kz2(-x');  similarly,  for  antisymmetric  mode 
we  letk;ji(x')= -kx2(-x')  and  k2i(x')=  -  k22('X')  since  the  currents  now 
have  an  odd  symmetry.  Then  the  coupled  I.E.  can  eventually  be  written  in 
terms  of  just  k;52(^')  and  k^(x').  We  will  solve  the  coupled  I.E.s  by  the 

method  of  moments  using  entire-domain  basis  functions. 


2-2  Entire-domain  basis  MoM  solution 

The  MoM  solution  to  equations  (10)  and  (11)  depends  on  a  reasonable 
choice  of  basis  functions.  These  basis  functions  should  incorporate  the  edge 

conditions  that  kz  is  singular  and  k^  is  zero  at  the  edges  of  the  microstrips.  In 
this  paper,  instead  of  using  subdomain  basis  MoM  for  which  the  solutions 
satisfy  the  electromagnetic  boundary  conditions  only  at  discrete  points,  we 
use  entire-domain  bases  for  expansion  of  the  unknown  surface  currents  and 
testing.  The  main  advantage  of  entire-domain  bases  lies  in  problems  where 
the  unknown  function  is  assumed  a  priori  to  follow  a  known  pattern  [9].  In 
such  cases,  entire-domain  functions  may  render  an  acceptable  representation 
of  the  unknown  while  using  far  fewer  terms  than  would  be  necessary  for 
subdomain  bases. 

Based  on  the  results  presented  in  [7],  it  is  found  that  Chebyshev 
polynomials  multiplied  by  an  appropriate  edge-factor  are  one  of  the  best  sets 
that  fit  the  solution  currents  of  this  coupled  microstrip  structure  well.  In  our 
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[ 

work  we  choose  Chebyshev  polynomials  of  the  second  kind  Un  (x')  to  expand 
the  unknown  transverse  transmission  line  currents  since  the  latter  can 
accurately  be  represented  with  a  linear  combination  of  just  few  functions  of 
this  kind.  By  similar  reason,  the  same  polynomials  of  the  first  kind  Tn(x') 
have  been  used  for  axial  currents  expansion.  Use  of  these  functions  has  the 
added  advantage  of  allowing  spatial  integrals  to  be  evaluated  in  closed  form. 


From  above  discussion,  we  expand  k^iC  x’)  and  kz2(  x’)  as  follows: 


k,2(x 


V  _  tT  /X'-b.^ 
n  =  0 


N-I 


kj2(x’)s 


V 


n  =  0 


where  an  and  bn  are  expansion  coefficients  representing  the  contribution  of 
each  order  of  the  Chebyshev  polynomials  U„(x')  and  T„(x')  to  the  unknown 
surface  currents  and  the  factors  in  front  of  the  summation  terms  give  U„(x’) 
and  Tn  (x')  the  correct  edge  behavior. 

To  test  the  coupled  I.E.s  after  doing  the  expansion  of  unknown  currents, 
we  follow  Galerkin's  method  and  utilize  the  same  functions  for  testing.  We 
need  the  tabulated  integral  formulas  listed  in  [10]  and  in  the  appendix  for  the 
calculation  involved  in  current  expansion  and  testing.  For  the  sake  of 
brevity,  we  omit  some  purely  algebraic  steps  and  finally  obtain  the  I.E.’s  for 
the  coupled  stmcture  as 
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symmetric  modes: 


^  N-l 

x:  (n+1) 

11=0 


^ —  cos  (4b)  Re  (j")  -  sin(4b)  Im  (j") 
^  Pc 


k;R-4  (R-c) 


Jm+l(^w)  Jn^i(4w)  d4 


N-l 

-Zb 

n=0 


n W  C  ^  ^  ^  ^m+l(^w)  J„(4w)  d4  =  0 


i 


'  [cos(5b)Re{j"}-sin(^b)Im{j”)](R-C')I„fiw)  J„i(|w)d4 


z:  5;  a„  (ii+l)C 

n=0  I  Pc 


ZbnW  I  •^— [cos(4b)Re{j"}-sin(4b)Im{j"}]  kcR-C  (R-C') 
n=0  I  Pc 


Jj4w)  J„(4w)d4=:0 


antisymmetric  modes: 


N-l 

z 

11=0 


8 

cr 

cos(4b)  Im  (j”)  +  sin  (4b)  Re 

k^R-4V-C') 

Jm+l(^w)  J„^i(4w)d4 


t>„  w  C  I  j^[cos(5b)  Im{j"}+siiifib)  Re{j"}](R-C)J„,,(5w)J„(4w)d^=:0 
n=0  J  Pc 
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The  integration  limits  (  —  00,00)  in  the  variable  ^  can  be  reduced  to 

[0, 00)  by  considering  m  and  n  to  be  even  or  odd  respectively.  Those  terms 

which  are  odd  with  respect  to  ^  will  integrate  to  zero  ,  and  those  even  terms 
will  survive  and  the  final  LE.s  can  be  written  as  two  times  the  same  integral 

with  the  integration  limits  (0,oo).  In  matrix  form  the  LE.s  can  be  written 
as 


2Nx2N 


ao 

• 

• 

• 

• 

^N-l 

• 

0 

bo 

• 

0 

• 

.  0, 

2Nxl  2Nxl 
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where  each  element  in  the  Z  matrix  is  in  terms  of  a  spectral  integral  on  and 

ao—.a^.  1,  bo,...,bN.  i  are  coefficients  for  expansion  polynomials.  We  have 
2N  equations  for  2N  unknowns.  To  obtain  a  non-trivial  solution,  the 
determinant  of  the  coefficient  matrix  must  vanish.  Since  Matrix  elements 
depend  on  the  propagation  constant  we  can  iterate  to  find  ^^r  the 
eigenmodes.  Some  typical  calculated  results  will  be  shown  in  next  section. 
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3.  Numerical  results 


A  computer  program  in  Fortran  77  using  the  method  described  above  was 
prepared  to  calculate  the  propagation  constants  (or  eigenvalues)  and  currents 
for  the  fundamental  and  higher  order  modes  of  the  coupled  microstrip 
structure.  For  convenience  in  comparison  with  the  results  in  [14]  and  [15], 
the  case  with  a  dielectric  constant  of  nf=kf/ko=3.13,  n(.=  1.0, 

t  =  0,635  mm,w  =  1.5  mm  and  a  varying  distance  b  has  been  calculated  and 
will  be  presented  here. 


3-1  Eigenvalues 

Table  1  shows  the  eigenvalues  for  the  fundamental  and  first  two  higher 
order  modes  at  operating  frequencies  near  cutoff.  For  each  mode  at  the 
specific  frequency,  it  can  be  seen  that  diey  agree  with  the  well  known  fact  that 

where  ^a,  ^iso  and  are  eigenvalues  corresponding  to 
antisymmetric,  isolated,  and  symmetric  modes,  respectively.  The  smaller  the 
spacing  b/w,  the  larger  (smaller)  the  eigenvalues  for  symmetric 
(antisymmetric)  system  modes.  When  the  separation  of  the  two  coupled 
microstrips  is  large  enough  (in  our  work,  about  b/w  >  2),  the  eigenvalues 
converge  to  the  corresponding  isolated  one.  Indeed,  the  symmetric  and 
antisymmetric  modes  can  be  viewed  as  emerging  from  the  corresponding 
isolated  mode.  Figure  3  shows  the  case  of  EH^  Mode.  While  (b-w)/w  is 
getting  smaller,  eigenvalues  of  the  EHq  symmetric  and  antisymmetric  system 
modes  are  separated  further,  and  ultimately,  when  two  microstrips  contact 
on  inside  edges,  the  EHq  symmetric  mode  goes  to  the  EH^  mode  of  an  isolated 
strip  with  double  width,  and  the  EHq  antisymmetric  mode  goes  to  the  EHj 
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mode  of  the  same  width-doubled  strip. 


The  behavior  of  eigenvalues  for  this  coupled  structure  agrees  very  well 
with  the  results  obtained  by  using  the  analysis  in  [1]. 


3-2.  Currents 

The  axial  and  transverse  currents  for  each  mode  mentioned  above  are 
shown  in  figures  4,  5  and  6  (since  currents  have  either  even  or  odd 
symmetry,  only  and  for  each  mode  are  shown  here.).  In  numerically 
quantifying  the  currents,  we  found  that  three  Chebyshev  polynomials  are 
enough  to  accurately  represent  the  unknown  surface  currents  on  an  isolated 
microstrip.  However,  for  currents  on  each  microstrip  of  the  coupled  system, 
we  need  the  first  five  to  obtain  an  accurate  result.  It  can  be  seen  in  Figure  4- 
a,  b,  5-b  and  6-b  that  due  to  the  repulsion  of  the  surface  charge,  currents  for 
symmetric  modes  have  smaller  value  near  the  inside  edge  than  at  the  outside 
edge.  On  the  other  hand,  due  to  the  attraction  of  charge,  the  antisymmetric 
modes  have  larger  current  near  the  inside  edge.  This  behavior  of  currents  is 
not  so  obvious  in  figure  5-a  and  figure  6-a,  since  at  such  low  frequencies,  the 
effective  width  of  the  strip  is  so  narrow  that  it  can  not  fit  the  complete  pattern 
of  the  currents.  At  higher  frequency  (40  GHz  for  EHj  mode  or  90  GHz  for 
EH2  mode,  for  example)  we  did  find  that  transverse  currents  for  EH^  a.nd 
EH2  have  the  same  behavior  as  mentioned  above.  In  all  cases,  we  can  see  the 
currents  for  isolated  line  lie  between  the  two  groups  of  currents  of 
symmetric  and  antisymmetric  system  modes  as  expected. 
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Conclusions 


A  full-wave  spectral-domain  integral  equation  formulation  has  been 
presented  for  analysis  of  coupled  microstrip  transmission  lines.  A  method  of 
moments  solution  utilizing  entire-domain  basis  functions  which  incorporate 
appropriate  edge  conditions  for  transverse  and  longitudinal  current 
components  allows  for  closed-form  evaluation  of  spatial  integrals.  In 
contrast  with  earlier  subdomain  basis  solutions,  improved  accuracy  is 
obtained  using  far  fewer  terms.  Numerical  results  in  the  form  of 
propagation  constants  and  current  distributions  for  the  dominant  and  first 
two  higher-order  modes  compare  favorably  to  results  of  other  techniques. 

A  more  difficult  problem  is  the  solution  in  the  case  of  lossy  modes.  This 
condition  occurs  when  the  propagation  constant  falls  below  the  eigenvalue  of 
a  surface-wave  mode  supported  by  the  integrated  circuit  background 
slabguide  structure.  In  this  case  singularties  in  the  spectral  integrands  move 
close  to  the  real  axis,  and  complex-plane  evaluation  techniques  become 
necessary  .  Work  along  these  lines  is  proceeding,  and  is  expected  to  yield 
results  in  the  near  future. 
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Appendix 


Other  than  the  tabulated  integral  formulas  listed  in  [10],  we  need  to 
develop  formulas  for  evaluating  following  integrals: 


cosax  U„(x)dx 


...  (1) 


sin  ax  Un(x)dx 


(2) 


1 


cosax  T^Cx) 


..  (3) 


To  evaluate  (1)  and  (2),  refer  to  Gradshteyn's  book,  eq.  (7-321),  p.  830, 
which  is 


(1-xV 


1 

2  iax 
e 


C/(x)dx  = 


r(2v+n) 
n!  r(v) 


a  3v+n(a) 


Since  U„(x)  =  C^i^(x),  where  C„(x)  is  the  Gegenbauer  polynomial.  The  above 
formula  can  be  specialized  to 


iax 

e  U„(x)dx 


K  i  r(n+2) 
n!r(l) 


-1 


a 


n 

i 


7c(n+l) 

a 
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Thus  we  obtain  following  formulas  for  (1)  and  (2); 


From  this  we  can  write 


for  odd  n. 


for  even  n. 
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1 


sinax  Tn(x) 


dx  _ 


{ 


(-1)2  JcJ„(a). 


0, 


for  odd  n . 

for  even  n. 


-1 


Also  notice  that 


n-> 

i  (-1)  2  7cJn(a),  for  odd  n. 

(-1)"^*7C  J„(a) ,  for  even  n. 


=  i  7cJ„(a), 


for  any  n . 


Thus  we  obtain  folowing  integral  fonnulas  for  use. 


/' 

*'-1 


cosaxTn(x)  ;  ^  =  Re  {  i  JcJn(a)  }  =7cJn(a)  Rc(i  } 

1 1-  x^ 
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1 

sinax  T„(x) -=^==Im{  i  tc  J„(a)  ]  =  tc  J„(a)  Im{i  } 
Y 1-  x-^ 

1 
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Table  1:  Eigenvalues 


^'^V.^Mode 

10  GHz 
EHq 

20  GHz 
EH^ 

40  GHz 
EH2 

1.10 

3.00638 

2.17164 

2.23297 

1.25 

2.98567 

2.12381 

2.21197 

1.50 

2.95200 

2.09482 

2.19654 

Isolated 

2.89586 

2.06867 

2.17533 

1.50 

2.83309 

2.05419 

2.15925 

1.25 

2.78686 

2.01026 

2.11995 

1.10 

2.75250 

1.94495 

2.05663 

I  Symmetric 
1  modes 


L  Antisymmetric 
modes 
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Figure  captions  : 


1.  Figure  1:  General  system  of  N  coupled  microstrip  lines. 

2.  Figure  2:  Two  identical  thin  coupled  microstrip  lines. 

3.  Figure  3:  Dependence  of  propagation  phase  constant  for  EHq  system 

modes  of  two  coupled  microstrip  lines  upon  line  spacing; 
symmetric  and  antisymmetric  system  modes  emerge  from  the 
isolated  EHq  mode. 

4.  Figure  4:  (a)  Transverse  currents;  (b)  axial  currents  for  EHq  coupled 

mode. 

5.  Figure  5:  (a)  Transverse  currents;  (b)  axial  currents  for  EH^  coupled 

mode. 

6.  Figure  6:  (a)  Transverse  currents;  (b)  axial  currents  for  EH2  coupled 

mode. 


cover  medium,  p 


oo 


Fig.l 


cover  medium, 


oo 


Normalized  Phase  Constant  (/k 


0.00  1.00  2.00  3.00  4.00 

Distance  between  coupled  strips  (b-w)/w 

Figure  3  ,  Dependence  of  propagation  _  phase  constant  (  for  EHo  system  modes 
of  two  coupled  microstrip  lines  upon  line  spacing;  symmetric  and 
antisymmetric  system  modes  emerge  from  the  isolated  EHq  mode. 
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Normalized  current  A^.n-l  Normalized  current  k^/k 


TRANSVERSE  CURRENTS.  COULPED  MODE  EHq 


LONGITUDINAL  CURRENTS.  COUPLED  MODE  EHq 


(n,-3.13.  nj-1.  t/w-.423.  f-10  GHz,  N-5) 


Position  X  (mm) 

Figure  4.  (b) 
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Normalized  current  A^.n-i  Normalized  current  k^/k 


TRANSVERSE  CURRENTS,  COUPLED  MODE  EHi 


(0) 


LONGITUDINAL  CURRENTS,  COUPLED  MODE  EH^ 


Position  X  (mm) 

Figure  5.  (b) 
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Normalized  current  A^^n-i  Normalized  current 


C  0.3 


TRANSVERSE  CURRENTS.  COUPLED  MODE  EH2 

(n,-3.13.  rtg-l.  t/w-.423.  f-40  GHz.  N-6) 


0:  Isolated  EH, 

\  symmetric  antisymmetric  b/w 


I  ■  '  ■  '  i  '  '  '  ‘  I  '  '  '  '  1 

0  .25  .5  .75  1.0  1.25  1.5  1.75  2.0  2.25  2.5  2.75  3.0 


Position  X  (mm) 

(a) 


LONGITUDINAL  CURRENTS,  COUPLED  MODE  EH2 

(nf-3.13,  Oe-I,  t/w-.423,  f-40  GHz,  N-6) 

-  \  symmetric  antisymmetric  b/w 


0:  isolated  EH. 


Figure  6. 
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Position  X  (mm) 
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4.6  EM  response  of  microstrip  patch  devices  and  antennas 

A  rigorous  description  for  natural  resonance  modes  of  microstrip  patch 
devices  in  the  MMIC  environment,  and  the  excitation  of  those  modes,  has  been 
completed  [16]  based  upon  an  EFIE  description  of  such  devices  immersed  in  a 
layered  dielectric/ferrimagnetic  surround.  Natural  resonances  arise  from 
simple-pole  singularities  in  the  temporal  complex  transform  plane,  and  are 
quantified  as  the  non-trlvial  solutions  to  resulting  homogeneous  integral- 
operator  equations.  Numerical  quantification  of  those  modes  for  a  circular- 
patch  device  is  implemented  [6]  using  Chebyshev  series  representations,  appro¬ 
priately  modulated  by  square-root  edge  singularity  factors,  for  the  induced 
current  components.  Complex  natural  frequencies  of  several  resonance  modes 
and  the  associated  natural-mode  currents  are  obtained.  A  M. S.  degree  thesis 
on  this  topic  is  presently  in  preparation. 

An  EFIE-based  coupled-mode  perturbation  theory  for  the  system-mode  reson¬ 
ant  frequencies  and  relative  current  amplitudes  describing  an  ensemble  of 
coupled  microstrip  antennas  in  the  MMIC  environment  [16]  has  been  generalized 
to  permit  the  presence  of  ferrimagnetic  film  layers.  This  perturbation  ap¬ 
proximation  exploits  natural-mode  currents  of  isolated  devices  to  obtain  coup¬ 
ling  coefficients  for  system-mode  currents  as  overlap  integrals  of  the  current 
on  one  device  with  the  resonant  field  of  another. 

Scientific  paper  [16]  was  included  with  Section  4.3,  while  paper  [6]  is 
appended  here. 
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Research  Paper  No.  [6] 


INTEGRAL  EQUATION  FORMULATION  FOR  NATURAL  MODES  OF  A  CIRCULAR 
PATCH  ANTENNA  IN  A  LAYERED  ENVIRONMENT 

Eric  W.  Blumbergs,  Dennis  P.  Nyquist,  and  Paul  F.  Havala 

Department  of  Electrical  Engineering  and  Systems  Science 
Michigan  State  University,  East  Lansing,  Michigan  48824 


This  paper  investigates  the  natural  resonant  modes  of  a  circu¬ 
lar  path  device  immersed  in  an  integrated  conductor,  film,  and 
cover  environment.  A  full  two-dimensional  EFIE  description  which 
totally  accounts  for  the  layered  environment  is  developed.  This 
development  uses  a  polar  coordinate  Sommerfeld- integral  representa¬ 
tion  of  the  dyadic  Green's  Function.  The  possibility  of  surface 
waves  of  the  layered  background  is  also  included. 

The  natural  mode  formulation  is  based  on  a  pair  of  coupled 
EFIE's.  These  equations  are  of  the  following  form: 

fJ=-ao  V=P,0  0 

for  asp, a  where  a^  (p ' )  are  radially  dependent  Fourier  expansion 
coefficients  of  surface  currents  on  the  patch,  K^^yCpPp')  are. ker¬ 
nels  arising  from  the  Hertzian  potential  Green's  dyad,  and  E^(p,0) 
are  components  of  the  impressed  field  tangent  to  the  patch  surface. 
Natural  modes  satisfy  the  homogeneous  specialization  of  eJ(p,0)=O. 

Orthogonality  of  the  exp(jn9)  leads  to  independent  systems  of 
homogeneous  equations  for  each  n.  Numerical  solutions  for  natural 
frequencies  and  eigenfields  are  pursued  by  Galerkin's  method.  The 
special  case  n=0  leads  to  independent  IE's  for  a  and  aig;  corres¬ 
ponding  field  components  are  coupled  for  higher-order  modes.  Rep¬ 
resentative  results  are  presented. 
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4.7  Experimental  measurement  of  EH  interactions  in  MMIC’s 

A  complete  V-band  mm-wave  instrumentation  setup  has  been  acquired  and 
installed.  Funds  provided  in  the  10/1/87  to  12/31/88  reneral  period  allowed 
the  purchase  of  a  required  local  oscillator  for  the  harmonic-mixer  detector 
and  MS-DOS  PCB  design  software  for  use  on  the  personal  computer  in  MSU’ s  new 
millimeter-wave  laboratory.  Low-loss  PCB  microstrip  circuits  are  being  fabri¬ 
cated  to  simulate  the  MMIC  environment.  Interactions  between  microstrip  waves 
and  adjacent  patch  devices  are  being  studied  experimentally  to  quantify  the 
associated  coupling  and  radiative  scattering.  The  measured  results  provide 
for  validation  of  the  analytically  predicted  scattering  characteristics.  Val¬ 
idation  of  the  coupled-mode  perturbation  theory  for  adjacent  microstrip  lines 
will  be  accomplished  through  experimental  measurements  on  PCB  implementations 
of  such  a  coupled-line  configuration. 
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